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Preiafig 

The  purpose  of  this  study  was  to  develop  and  to  demonstrate  a 
'  control  theory,  incorporating  the  feedback  of  fractional  order 
derivatives  of  the  structural  response,  for  the  control  of  undamped  and 
viscoelastically  damped  structures.  The  control  theory  was  developed 
from  the  finite  element  formulation  of  the  structure,  and  the  use  of 
fractional  order  state  feedback  provided  control  authority  equal  to 
current  structural  control  schemes  with  the  use  of  fewer  sensors. 
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Abstract 

'U/A-iC 

The  purpose  of  this  study  is  to  demonstrate  the  analysis  leading  to 
the  development  of  the  fractional  order  state  feedback  control  theory 
for  structural  control  of  both  undamped  and  viscoelastically-damped 
structures.  It  is  shown  herein  that  there  exists  a  relation  between  the 
traditional  state  vector  that  includes  structural  displacements  and 
velocitie8f  and  the  fractional  state  vector  which  includes  fractional 
derivatives  of  structural  responses.  This  relation  permits  the 
modification  of  linear  quadratic  regulator  theory  to  include  the 
application  of  fractional  order  states  in  the  feedback  control.  The 
application  of  this  theory  leads  to  an  alternative  form  of  an  observer. 


FRACTIONAL  ORDER  STATE  FEEDBACK  CONTROL  FOR  UNDAMPED  AND 


VISCOELASTICALLY-DAMPED  STRUCTURES 


I.  Introduction 


Structural  control  of  large  space  structures  has  become  of  interest 
in  recent  years  due  to  the  application  of  these  structures  in  areas  such 
as  the  NASA  space  station  and  the  Strategic  Defense  Initiative.  The 
large  size  and  low  weight  of  these  structures  result  in  large  number  of 
closely  spaced,  low  frequency  vibration  nodes,  which  require  structural 
control  to  aaintain  systen  stability  and  performance  requirements. 
Passive  damping  systems  such  as  damping  pads  and  sophisticated  shock 
absorbers  provide  only  a  partial  solution  to  the  problem  due  to  their 
limited  control  action  at  low  frequencies.  The  incorporation  of  active 
structural  control  systems  with  passive  damping  devices  shows  potential 
in  improving  overall  systen  performance  considerably. 

The  current  research  effort  began  with  an  interest  to  introduce 
active  control  to  structures  that  are  viscoelastically  damped,  in 
searching  for  improved  system  response.  The  objective  of  this  thesis  is 
to  demonstrate  the  analysis  leading  to  the  development  of  a  fractional 
order  state  feedback  control  law,  for  structures  incorporating  both 
passive  damping  through  viscoelastic  materials  which  are  modeled  by  a 
fractional  derivative  stress-strain  relation  and  active  damping  by 
applied  control  forces.  To  achieve  this,  quadratic  optimal  control 
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theory  is  aodified.  Specifically,  linear  quadratic  regulator  theory  is 
modified  to  include  fractional  derivatives  in  the  state  vector.  It  is 
shown  herein  that  there  exists  a  relation  between  the  traditional  state 
vector  and  the  fractional  order  state  vector  which  permits  the 
incorporation  of  fractional  order  states  in  the  feedback  control  law. 

The  thesis  first  reviews  the  properties  of  the  quadratic  optisal 
control  theory  and  the  fractional  derivative  in  the  formulation  of  a 
viscoelastic  constitutive  law.  The  development  of  the  fractional  order 
state  feedback  control  theory  is  presented  in  Chapter  III  for  both 
undamped  and  viscoelastically-damped  structures.  Two  example  problems 
in  Chapter  IV  illustrate  the  application  of  this  control  theory  and 
demonstrate  the  solution  technique. 


II.  Brief  Overviews  of  Optimal  Control  Theory  and  Fractional  Derivatives 


Before  proceeding  with  the  development  of  the  fractional  order 
state  feedback  control  theory,  it  is  appropriate  to  introduce  the 
properties  of  the  optimal  control  theory  and  fractional  derivatives 
relevant  to  the  following  control  theory  development.  Of  particular 
interest  are  the  linear  quadratic  regulator  and  the  basic  fractional 
derivative  viscoelastic  constitutive  relation.  Section  2.1  provides  a 
brief  overview  on  the  development  of  the  linear  quadratic  regular 
control  theory.  Section  2.2  provides  a  general  description  of  the 
fractional  calculus  and  the  development  of  the  fractional  derivative 
constitutive  relation  for  viscoelastic  materials. 

2.1  Optimal  Control  Systems  Based  on  Quadratic  Performance  Index 

The  concept  of  control  system  optimization  compromises  the 
selection  of  a  performance  index  and  a  control  law  which  yields  the 
optimal  control  system  within  the  limits  imposed  by  system  physical 
constraints.  A  performance  index  is  a  function  whose  values  indicates 
how  well  the  actual  system  performance  matches  the  desired  system 
performance.  In  many  cases,  system  behavior  is  optimized  by  choosing 
the  control  vector  in  such  a  way  that  the  performance  index  is  minimized 
(or  maximized). 

In  many  practical  control  systems,  it  is  desired  to  minimize  a 
measure  of  the  difference  between  the  actual  state  and  the  desired 
state,  and  to  minimize  the  energy  required  for  the  control  action.  For 
a  given  linear  system  represented  by 


x  =  Ai  +  Bu 


(1) 


where 

£  =  state  vector  (n-diaensional  real  vector) 
u  a  control  vector  ( r-diaensional  real  vector) 

A  =  n  x  n  aatrix 
B  s  n  x  r  aatrix 

the  associated  quadratic  perforaance  index  of  the  systea  over  the  tiae 
interval  0  *  t  *  t  Bay  be  written  as  (1:753) 

>  T  Vt 

J  *  J  [C(t)  -  x( t ) ]  Q[C(t)  -  x( t ) ]  dt  +  xj  u  (t)Bu(t )dt  (2) 

o  0 


where 


C(t)  =  the  desired  state 
x(t)  a  the  actual  state 
u(t)  a  the  control  vector 

Q  =  a  positive-definite  (or  a  positive-aeaidefinite)  aatrix 
B  *  a  positive-definite  aatrix 

The  Lagrange  aultiplier,  X,  is  a  positive  constant  indicating  the  weight 
of  control  energy  with  respect  to  the  ainiaizing  errors,  and  the  tiae 
interval  0  £  t  £  t  is  either  finite  or  infinite. 
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Finding  the  optiaal  control  law  for  tne  systea  subject  to  the 
perforaance  index  given  by  Eq  (2)  has  a  practical  significance  that  the 
resulting  systea  coaproaises  between  ainiaizing  the  integral  error  and 
ainiaizing  the  control  energy.  For  the  case  where  tf=  and  the 
desired  final  state  is  the  origin,  or  C  s  0,  the  quadratic  perforaance 
index  given  by  Eq  (2)  can  be  expressed  as 

00 

J  =  J[x(t)TQx(t)  +  uT(t)Bu{t)]dt  (3) 

o 

where  X  has  been  included  in  the  positive  definite  aatrix  8.  This  type 
of  control  problea  is  often  referred  to  as  a  "linear  quadratic  regular 
problea”  since  the  final  state  response  is  desired  to  be  zero  (1:126). 

In  designing  optiaal  control  systeas  based  on  the  quadratic 
perforaance  index  shown  in  Eq  (3),  one  is  interested  in  choosing  a 
control  vector,  u(t),  so  that  the  perforaance  index  is  ainiaized.  The 
necessary  conditions  for  an  optiaal  solution  of  the  regulator  problea 
aay  be  obtained  by  using  the  calculus  of  variations  which  nay  involve 
the  concept  of  a  Haailtonian  (2:180-181).  Once  the  initial  and  final 
state  conditions  have  been  identified  by  the  state  and  costate 
equations,  the  optiaal  control  law  aay  be  deterained  by  solving  the 
two-point  boundary-value-problen  (2:189-190).  The  optiaal  control  for 
the  regulator  problea  is 

11  =  -  B"lBSx  (4) 

where  S  is  the  solution  of  the  algebraic  aatrix  Riccati  equation: 


(5) 


0  =  ATS  +  SA  -  SBB'Vs  +  Q 
Defining  the  optiaal  gain  aatrix  as 

G  =  R"Vs 

the  optiaal  control  can  be  written  as 


(6) 


u  =  -  Gx 


(7) 


and  the  state  equation  can  be  expressed  as 

k  =  Ag  -  BGx  (8) 

2.2  Fractional  Derivatives  and  Its  Application  to  Viscoelastic  Materials 

A  fractional  calculus ,  in  which  fractional  order  derivatives  and 
integrals  are  defined  and  studied,  is  not  a  very  new  concept.  Although 
the  idea  of  fractional  calculus  has  been  in  existence  since  the  late 
1600’s,  the  use  of  fractional  calculus  has  been  few  as  compared  to  the 
classical  calculus  of  integer  orders.  An  article  written  by  Ross  (3) 
provides  an  interesting  overview  on  the  historical  developsent  of  the 
basic  concepts  of  fractional  calculus. 

A  fractional  derivative  is  a  linear  operation  that  generalizes  the 
order  of  differentiation  to  fractional  values.  It  is  the  inverse 
operation  of  fractional  integration  attributed  to  Riesann  and  Liouville 
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(4).  The  fractional  derivative  is  defined  as 


Da[x(t)] 


1 

ra-a) 


x(T) 


( t“T  y 


dr 


(9) 


for  0  <  a  <  1.  The  fractional  derivative  operator  has  the  property  in 
the  Laplace  transform  doaain  that 


Sa-C[x(t)] 


(10) 


where 


00 

£[x(t)l  s  J  x(t)  e~at  dt  (ID 

o 

A  similar  relationship  exists  for  the  fractional  derivative 
operator  in  the  Fourier  transfora  doaain 

sjy'fxU)]]  *  (iw)a-3[x(t)]  (12) 

where 


CO 

*5[x( t )  1  *  J  *(t)  e  *Wt  dt  (13) 

o 

In  fact,  the  fractional  derivative  satisfies  many  of  the  saae 
properties  as  the  integer  derivative,  particularly  linearity  and  the 
composition  property 
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(14) 


Da[ay(t)  +  bx(t)]  *  aDa[y(t)]  *  bDa[x(t)] 

for  a  and  b  are  constants,  and 

Da[lftx<t)]]  =  Da+P[x(t)]  (15) 

The  development  of  the  fractional  derivative  stress-strain 
constitutive  model  for  viscoelastic  materials  has  its  foundation  in  the 
early  observations  made  by  Nutting  and  Geaant  (6,7,8).  Nutting's 
observations  suggested  that  the  stress  relaxation  property  could  be 
modeled  by  fractional  powers  of  time,  while  Gemant's  observations 
suggested  that  the  frequency  dependent  properties  of  viscoelastic 
materials  could  also  be  modeled  by  fractional  differentials  of  time. 

The  simultaneous  modeling  of  both  relaxation  and  frequency  dependent 
characteristics  of  viscoelastic  materials,  using  fractional  time 
differentials,  was  later  proposed  by  Scott-Blair  in  1950’s  (9). 
Meanwhile,  the  use  of  fractional  derivatives  was  suggested  by  Caputo  for 
modeling  the  viscoelastic  behavior  of  geological  strata  (10,11).  He  and 
Minard  also  showed  through  experimental  data  that  the  fractional 
derivatives  constitutive  relationship  could  be  used  for  some  glasses  and 
metals  (12).  Bagley  proposed  the  use  of  fractional  derivatives  in  the 
finite  element  analysis  of  viscoelastically  damped  structures  (5).  In 
recent  years,  Bagley  and  Torvik  have  jointly  published  several  papers 
demonstrating  the  application  of  fractional  derivatives  in 
viscoelasticity  and  structural  dynamics  (13-16).  They  have  also  shown 
that  the  fractional  order  models  are  consistent  with  thermodynamic 
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constraints  (17)  and  these  aodels  have  their  foundation  in  aolecular 
theories  that  predict  the  aacroscopic  behavior  of  viscoelastic  aaterials 
(18). 

The  general  fora  of  the  fractional  derivative  viscoelastic  aodel 
put  forward  by  Bagley  and  Torvik  is  (14) 

M  _  N 

<y(t)  +  2  b  DPV(t)l  =  Ee(t)  +  £  E  Dan[e(t)l  (16) 

hi  o  n 

r»  »  1  n  « 1 

where  the  tine-dependent  stress  fields  are  related  to  the  tiae-dependent 
strain  fields  through  series  of  fractional  derivatives.  For  aany 
viscoelastic  aaterials,  the  stress-strain  relation  can  be  aodeled  by 
retaining  only  the  first  fractional  derivative  tern  in  each  series  in  Eq 
(16)  and  the  resulting  aodel  is 

<H t)  +  bD^[<r(t)J  *  E  e(t)  +  E,Da[e(t)]  (17) 

O  1 

The  Second  Law  of  Theraodynaaica  requires  that  (17) 

E  i  0  E.  i  bE 

o  1  o 

Ej  >  0  a  *  0  (18) 

b  >  0 

These  constraints  ensure  that  energy  dissipation  rate  and  stored  strain 
energy  are  positive. 

The  stress-strain  relation  expressed  in  the  Laplace  doaain  is 


(19) 


Q(  s) 
E(8) 


G 

o 


,  ,  _a 

1  +  bS 


The  result  is  a  viscoelastic  sodel  with  four  parameters:  b,  E  ,  E<(  and 
at.  This  oodel  has  been  shown  to  be  very  accurate  over  several  decades 
of  frequency  for  approximately  130  saterials  (15,16).  For  sose 
applications  where  the  required  frequency  range  is  lisited,  a  three 
paraaeter  viscoelastic  model,  where  b  =  0,  is  sufficient  (21:129).  The 
three  parameter  model  is  defined  as 

■Qsl  s  E  +  E.S®  (20) 

e(s)  o  l 

and  this  is  the  model  that  will  be  used  in  this  work. 
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III.  Fractional  Order  State  Feedback  Control  Theory 


The  development  of  the  fractional  order  state  feedback  control 
stems  from  the  interest  in  examining  the  effect  of  feedback  in 
controlling  viscoelastically  damped  structures.  Bagley  and  Calico  have 
shown  that  fractional  order  state  equations  can  be  developed  where  the 
state  vector  includes  fractional  order  time  derivatives  of  structure 
motion  (19).  Their  work  also  suggested  possible  improvement  in  system 
performance  by  feeding  back  fractional  time  derivatives  of  system 
response.  The  fractional  order  state  equations  for  a  general  system 
prescribed  by  Bagley  and  Calico  is  (19:493) 

D°z  =  Ay  -  BGy  (21) 

for  z  is  the  fractional  order  state  vector  and  0  is  the  matrix  of 
effective  gain  coefficients. 

In  order  to  apply  fractional  order  derivatives  feedback  in 
classical  structural  control,  it  is  necessary  to  develop  a  control  law 
incorporating  the  fractional  order  states.  The  development  of  the 
theory  will  begin  with  the  formulation  of  the  fractional  order  equations 
of  notion,  which  is  an  extension  of  the  traditional  equations  of  motion 
derive  from  the  finite  element  formulation  of  the  structure  of  interest. 
A  special  solution  technique  is  then  demonstrated  for  determining  the 
the  eigenstructure  of  the  fractional  order  equations  of  notion  which 
includes  both  integer  order  and  fractional  order  time  derivatives  of 
structural  displacements.  The  eigenvectors  of  the  system  are  then  used 
to  develop  the  special  traditional'-state/fractional-state  transformation 
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relation,  and  the  result  of  this  transformation  is  then  applied  to  the 
development  of  the  control  law  which  incorporates  the  feedback  of 
fractional  order  states. 

In  the  application  of  optimal  quadratic  control  theory  for  active 
structural  control,  more  specifically  for  the  linear  quadratic  regulator 
problem,  it  has  been  shown  that  the  quadratic  control  theory  produces 
optimal  control  as  linear  feedback  of  the  state  vector  (20:129).  The 
state  equations  for  the  structure  can  be  written  as 

S  s  As  -  BOj  (8) 

where  x  is  the  state  vector  that  includes  only  structural  displacements 
and  integer  order  tine  derivatives  of  structural  displacements,  and  Q  is 
the  optimal  gain  matrix.  It  will  be  shown  herein  that  there  exists  a 
relation  between  the  traditional  state  vector  of  integer  order  and  the 
fractional  order  state  vector  given  by 

X  ■  ♦  X  (22) 

where  ♦  is  as  a  state  transformation  matrix  which  can  be  determined 
from  the  eigenstructure  of  the  equations  of  motion  posed  in  an  expanded 
format. 

Consider  the  finite  element  equations  of  notion  for  a  structure  with 
N  degrees  of  freedom 

[N] (w(t) }  ♦  [K]{w(t)}  *  {P(t)>  (23) 
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where  [N]  and  [K]  are  the  lass  and  stiffness  aatrices  of  the  finite 
eleaent  aodel,  and  (F(t)}  is  the  forcing  function.  In  order  to  develop 
the  state  transforaation  aatrix,  <P,  it  is  necessary  to  first  pose  Eq 
(23)  as  a  set  of  fractional  order  differential  equations  and  then 
deteraine  the  eigenstructure  associated  with  the  fractional  order 
equations  of  aotion.  Applying  the  coaposition  property  of  the 
fractional  derivative  such  that 

D°t[l)P{w(t))]  *  Da+^{w(tU  (24) 

and  Eq  (23)  can  be  posed  as 

D®[M]jD2~®  +  D2-2®  +  D2*3®  +  •—  +  D°j  (w(t) ) 

-  [M]£d2"®  +  D2"2®  +  D2”3®  ♦  —  +  D®j (w(t) } 

+  [KjD°{w(t)>  =  (F(t) }  (25) 

Eq  (25)  represents  the  fractional  order  equations  of  aotion  for  the 
systeB  defined  in  Eq  (23).  Here  a  is  chosen  to  be  a  fraction  of  the 
fora  1/n  where  n  is  an  integer.  This  is  to  ensure  that  velocities  will 
appear  in  the  fractional  order  state  equations  and  thus  allows  one  to 
solve  the  initial  value  problea  where  the  initial  velocity  is  specified. 
For  systens  with  sero  initial  condition,  this  restriction  on  a  is  not 
required. 

The  eigenstructure  of  the  fractional  order  equations  of  notion,  Eq 
(25),  can  not  be  solved  using  traditional  eigenvalue  problea  solution 
aethod  because  of  the  inclusion  of  fractional  derivatives.  The  special 
solution  technique  presented  here  for  solving  Eq  (25)  is  an  extension  of 
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the  method  presented  by  Foss  for  determining  the  solution  of  an  N 
degree-of-freedom  (DOF)  viscously  damped  structure  (29:361) 

[M] (q(t) }  ♦  [C](q( t ) }  ♦  [K]{q(t) )  =  {F(t)>  (26) 

where  the  daaping  matrix,  [C],  is  not  a  linear  coabination  of  the  aass 
and  stiffness  matrices  of  the  system.  Foss  proposed  the  use  of  a  set  of 
auxiliary  variables  and  the  conversion  of  a  Lagrangian  set  of  N  second 
order  differential  equations  into  an  equivalent  set  of  2-N  first  order 
differential  equations  known  as  Hamilton’s  canonical  equations.  This 
leads  to  the  following  equivalent  equations  of  motion: 

[MHy(t)}  +  CK] ( y( t ) }  *  {Y(t)>  (27) 


where 


[Ml  * 

[0] 

[M] 

.  [Ml 

[Cl  . 

CK]  - 

-[Ml 

(01 

.  [01 

CK]  . 

(y(t)}  *  <  ► 

(q(t)> 


[Y( t) ) 


{‘°»  } 

L  {F(t)»  J 


(28) 


(29) 


(30) 


(31) 
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The  lower  half  of  the  partitioned  aatrix  equations,  Eq  (27),  is  the 
equations  of  aotion  of  the  systea  in  Eq  (26)  and  the  upper  half  of  the 
aatrix  equations  is  satisfied  identically.  The  foraulation  of  Eq  (27) 
has  the  advantage  that  its  solution  can  be  deterained  using  traditional 
eigenvalue  problea  solution  aethod. 

To  deteraine  the  eigenstructure  of  Eq  (25),  one  can  adopt  the 
aethod  presented  by  Bagley  and  pose  the  fractional  order  equations  of 
aotion  in  teras  of  two  real,  square,  syaaetric  aatrices  in  the 


following  foraat  (14:744) 

Da[Ml [w(t) )  +  [K] {w( t) }  ={F(t) } 

(32) 

CO] 

[01 

—  [01 

[Ml" 

[01 

[0] 

•—  [Ml 

[01 

CM]  * 

j 

1  | 

j 

(33) 

[01 

[Ml 

—  [0] 

[01 

[Ml 

[01 

—  [01 

[03 

[01 

[01  - 

-  [0] 

-[Ml 

[01 

[01 

[0J  - 

~  -[Ml 

[01 

[01 

[K]  *  i  i  I 


[01 

-[Ml 

[0] 

[01 

[01 

-[M] 

[01 

—  [01 

[01 

[01 

[01 

[01 

~~  [0] 

[01 

CK] 

(34) 
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'  D(2'a) 
~<2-2a) 


(w(t))  =  < 


.2a 


jx 


{w(t) }  ) 
{w(t> } 

{W( t>  > 

<w(t)>  J 


(  {0} 
(0) 


<F(t) >  =  { 


10} 

UF(t))J 


(35) 


(36) 


This  revised  for*  is  referred  to  *s  the  expanded  equations  of  sotion 
where  [MI  and  [K]  are  the  pseudo-sass  and  pseudo-stiffness  aatrices 
respectively  (15:922).  Note  that  the  lowest  set  of  partitioned  satrix 
equation  is  the  equations  of  notion  given  in  Eq  (25),  and  that  all  of 
the  upper  sets  are  satisfied  identically.  The  order  of  the  expanded 
equations  of  notion  is  2-N*n. 

Setting  (P(t)}  »  {0}  yields  the  honogeneous  forn  of  Eq  (32)  and  the 
expanded  eigenvalues  and  eigenvectors  for  the  systen  nay  now  be  solved 
to  satisfy  the  equation 


+  [KHfy  *  {0}  (37) 

The  solutions  to  this  equation  can  be  conputed  using  one  of  several 
techniques  currently  available  (22).  The  eigenvalues  appear  in  conplex 
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conjugate  pairs  and  there  are  Nn  pairs  of  eigenvalues  for  the  systea. 
The  eigenvalues  nay  be  aapped  onto  the  Rieaann  surface  for  the  function, 

^  1  /  fl 

Xj  =  Z  ,  consisting  of  n  Rieaann  sheets  (23:902).  N  pairs  of  coaplex 
conjugate  eigenvalues  usually  appear  to  each  sheet.  The  eigenvalues  on 
the  principal  sheet  of  the  Rieaann  surface,  which  is  the  first  upper 
sheet  of  the  n  sheets  present,  represent  poles  in  the  systea  transfer 
function  and  they  produce  sinusoidal  notion  of  the  structure  at  the 
resonant  frequencies  of  the  systea.  The  eigenvectors  associated  with 
the  eigenvalues  on  the  principal  sheet  describe  the  node  shapes  of  the 
structure  at  the  nodal  frequencies.  The  eigenvalues  on  other  sheets  of 
the  Rieaann  surface  represent  polea  in  the  systea  transfer  function 
which  produce  a  aonotonically  decreasing  response  of  the  structure. 

This  aonotonically  decreasing  notion  describes  the  creep  and  relaxation 
response  of  the  systea.  It  will  be  shown  in  Chapter  IV  of  this  thesis 
that  the  systea  response  is  strongly  influenced  by  the  residues  of  the 
poles  on  the  principal  sheet  of  the  Rieaann  surface,  and  one  aay 
describe  the  general  response  of  the  systea  with  very  good  accuracy  over 
an  internediate  range  of  tine  using  only  the  eigenvalues  and 
eigenvectors  associated  with  those  poles  on  the  principal  sheet.  In  the 
following  analysis,  only  the  eigenvalues  and  eigenvectors  associated 
with  the  poles  on  the  principal  sheet  of  the  Rieaann  surface  are 
retained  for  the  foraulation  of  the  state  transforaation  aatrix. 

Recalling  the  expansion  theoren  in  structure  nodal  analysis 
(24:283)  where  the  systea  response  aay  be  described  as  a  superposition 
of  the  noraal  nodes  of  the  systea  aultiplied  by  corresponding 
tine-dependent  generalized  coordinates  and  be  written  in  the  fora 


17 


{q}  *  [uj{n> 


(38) 


a  a  coluan  vector  of  the  generalized  systea  response 
a  the  aodal  aatrix  consisting  of  the  aodal  vectors  arranged  in 
a  aatrix 

a  a  coluan  vector  of  the  aodal  coordinates  of  the  systea 


Applying  the  theorea  to  the  expanded  systea,  one  aay  now  pose  the  systea 
response  of  the  structure  in  the  fora 


{W(t)>  a 


where 


D(Z_a>  {w(t)> 
D(2'ZO){w(t)} 


{w(t)>  a 


{w(tn 


{"UU 

{w(t) } 


[♦]  *  uv  cv  fi2)  {5mh 
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(41) 


Notice  the  expanded  response  (w(t)}  includes  nodal  displaceaents, 
velocities,  and  the  fractional  derivatives  of  the  nodal  displaceaents 
for  all  N  degrees  of  freedoa  of  the  finite  eleaent  aodel.  The  aodal 
aatrix,  [4],  consists  of  the  conjugate  pairs  of  aodal  vectors,  and 

{$!>,  associated  with  the  conjugate  eigenvalues,  X^  and  X^ ,  on  the 
principal  sheet  of  the  Rieaann  surface. 

Froa  Eq  (39)  one  nay  extract  the  necessary  inforaation  to  construct 
the  traditional  state  vector  and  the  fractional  state  vector  for  the 
systea.  The  traditional  state  vector  aay  be  posed  in  a  fora  similar  to 
Eq  (39) 

X  *  [*AHn}  («) 

where 


X 


'  {w(t)) 

k  D,{w(t)}  . 


(43) 
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[♦Al  a  a  aatrix  consists  of  row  vectors  in  Eq  (40)  that  are 
associated  with  the  states  defined  in  Eq  (43) 


Siailarly  the  fractional  state  vector  say  be  written  as 


x  =  [$BHn} 


where 


(44) 


(45) 


[♦Bl  «  a  aatrix  consists  of  row  vectors  in  Eq  (40)  that  are 
associated  with  the  states  defined  in  Eq  (45) 

The  fractional  state  vector  in  Eq  (45)  is  a  special  subset  of  the  state 
vector  of  the  expanded  equations  of  notion  in  Eq  (35)  and  each  state  of 
the  fractional  state  vector  represents  a  state  of  interest  in  Eq  (35). 
The  fractional  state  vector  in  Eq  (45)  has  the  saae  dinension  as  the 
traditional  state  vector  in  Eq  (43)  and  it  includes  a  total  of  2-N 
fractional  and  integer  states.  The  states  of  the  fractional  state 
vector  nay  be  of  any  integer  or  fractional  order  states  taken  froa  one 
structural  node  or  froa  a  conbination  of  several  different  nodes. 
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Fro*  Eq  (44),  one  may  define  {17}  as 

{r?}  *  [* BI_1z 

Substitute  Eq  (46)  into  Eq  (42)  yields 

X  * 

Define  the  state  transformation  satrix,  <t>,  as 

♦  *  rojEV'1 

and  Eq  (47)  can  be  written  as 

X  *  *  X 


(46) 


(47) 


(48) 


(49) 


Eq  (49)  illustrates  the  relation  between  the  traditional  state 
vector  and  the  fractional  state  vector  in  a  system,  and  the  state 
transformation  matrix  captures  the  essence  of  this  important  relation 
(19).  The  task  at  hand  is  to  utilize  the  traditional-state/fractional- 
state  relation  just  presented  as  the  foundation  for  developing 
fractional  order  state  feedback  control  law  for  both  undamped  and 
viscoelastically  damped  structures.  The  development  of  the  theory  will 
be  presented  in  the  following  sections. 

1 
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3.1  Fractional  Order  State  Feedback  Control  of  Undamped  Structures 


In  this  section,  the  fractional  order  state  feedback  control  law 
will  be  developed  for  the  control  of  undaaped  structures.  The  optimal 
quadratic  regulator  control  theory  serves  as  the  starting  point  for  the 
developaent  of  the  new  control  law.  The  traditional-state/fractional- 
state  relation  is  then  incorporated  into  the  theory  for  the  formulation 
of  the  factional  order  state  feedback  control  law. 

Consider  the  equations  of  motion  of  an  N  degree  of  freedoa 
structure  under  active  control 

[M]{w(t)}  ♦  [K]{w(t) }  *  [Dl{u(tH  (50) 

The  equations  of  aotion  for  the  system  are  similar  to  Eq  (23),  except 
for  the  forcing  function  which  ham  been  replaced  by  the  product  of  a 
matrix,  [0],  and  a  vector,  (u(t)}.  The  matrix, [D],  describes  the 
location  and  orientation  of  the  actuators,  and  the  vector,  (u(t)}, 
represents  the  actuator  force.  The  product  of  [0]  and  (u(t)}  represents 
the  control  force  applied  to  the  structure. 

The  system  in  Eq  (44)  may  be  expressed  in  state  model  fora 

i  *  A*  ♦  By  (1) 


where 


X  * 


(51) 
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(52) 


(53) 


(54) 


and  [0]  is  a  zero  aatrix  and  [I]  is  an  identity  matrix. 

Recalling  fros  optisal  quadratic  regulator  control  that  the  optisal 
for  the  systes  is  a  linear  feedback  of  the  state  vector 


U  =  -  Gx  (7) 

for  G  is  the  optisal  gain  matrix  for  the  system  derived  from  the 
solution  of  the  Riccati  equation.  Introducing  the  optimal  control  into 
the  system  and  the  state  equations  nay  be  written  as 

X  ~  Ax  -  BGx  (8) 

Partition  the  gain  matrix  G  into  two  equal  size  matrices,  Gt  and  G2, 
such  that 


G  *  [G,:G2J  (55) 

and  incorporate  the  gains  into  the  finite  element  model  of  the  structure 
to  obtain 
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(56) 


[M](w(t) }  +  CD] [02){w(t> >  +  UK]  +  [D][Gt]){w(t)>  *  (0) 

Note  that  the  feedback  control  introduces  control  daaping  and  control 
stiffness  to  the  structure,  as  evidenced  by  the  inclusion  of  [D][G2J  and 
[DUG,]  in  Eq  (56). 

To  determine  the  nodal  Matrix  and  the  state  transformation  aatrix 
for  the  structure  under  control,  first  pose  the  equations  of  notion,  Eq 
(56),  in  the  fractional  foraat 

D°‘tM][D2"a  +  D2‘20t  ♦  D2_30t  +  •—  ♦  D°|  (w(t) } 

-  (M] ^D2_a  ♦  D2-2a  ♦  D2_3a  ♦  --  ♦  D°| { w( t ) } 

+  d“[D][G2][d,-<X  +  D1"2®  +  D1"3**  +  +  D°j  (w(t) } 

-  [D](G2][d1'<X  +  D12<X  ♦  D1"3®  +  .  +  D®|  (w(  t) ) 

+  ([K]  +  [D][Gt])D°(w(t))  *  {0}  (57) 

The  eigenvalues  and  eigenvectors  of  the  systea  in  Eq  (57)  aay  now  be 
solved  using  the  solution  technique  of  the  expanded  equations  of  notion 
presented  in  Eq  (32)-Eq  (36).  Por  the  equations  of  notion  are  now 
written  as 

Da[M]{w(t))  +  [K] (w(t) )  *  { F( t ) )  (32) 
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where 


[0] 

[0] 

*—  "*** 

[01 

[Ml 

[0] 

[01 

— *  **— * 

— 

[Ml 

CM]  = 

!  j 

[dhg2i 

(58) 

(01 

[Ml 

— 

[dhg21 

[0] 

[01 

[Ml 

[01 

[DJ[G2J 

— 

[01 

[01 

’  [01 

[01 

..... 

— 

[0J 

-[Ml 

[0] 

[01 

[0] 

. 

— 

-tM] 

[0] 

[0] 

[K]  * 

j 

i 

2  2 

t  j 

j  - 

[D][G2] 

j 

(59) 

[01 

-[M] 

— 

-[d][g2] 

[0] 

[01 

[0] 

-[MJ 

[0] 

— [D][02 

1  -• 

[01 

[01 

[01 

[0J 

[0] 

[01 

[0] 

[KJ+[D][G  1 

and  { w( t ) }  is  the  state  vector  defined  in  Eq  (35)  and  (F(t)}  is  a  zero 
vector  for  this  case.  Following  the  procedure  identified  in  Eq  (39)-Eq 
(41),  the  aodal  aatrix,  [9],  say  be  detersined  for  the  expanded  response 
of  the  structure  under  feedback  control.  Furthermore  the  state 
transformation  aatrix  may  be  formulated  following  the  steps 
identified  in  Eq  (42)-  Eq  (48).  Applying  the  state  transformation 
relation, Eq  (49),  to  the  optimal  feedback  control,  Eq  (8),  and 

i  S  Ag  -  BG*X  (60) 

This  is  the  control  law  for  active  structural  control  incorporating 
the  feedback  of  fractional  order  derivatives  of  the  structure's 
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response.  The  use  of  fractional  state  feedback  control  peraits  one  to 
■aintain  equivalent  control  authority  over  the  systea  as  the  traditional 
state  feedback  control  with  the  use  of  fewer  sensors  (additional 
control  hardware  and  circuitries  say  be  required  to  support  the  actual 
operation).  The  additional  inforaation  provided  by  the  fractional 
derivatives  of  the  structural  response  and  the  use  of  the  state 
transforaation  aatrix  allows  one  to  observe  the  notion  at  all  structural 
nodes  of  interest,  even  though  only  United  sensors  are  actually  used  to 
aonitor  the  systea. 

To  further  expand  the  theory,  substitute  the  state  transforaation 
relation  into  Eq  (60)  and  the  control  nay  now  be  written  as 


<Pi  »  A<Pz  -  BG$y 


Multiply  both  sides  of  Eq  (61)  by  ♦_1  and 


(61) 


I  a  <P~1A<Pz  -  $"1BG$y 

Redefine  the  eleaents  of  Eq  (62)  as  follows 


(62) 


A*  =  ♦''A# 

_*  1_ 

B  a  <p  B 


G+ 


(63) 

(64) 

(65) 


and  Eq  (62)  aay  be  written  as 

I  *  A  i  -  0  i 


(66) 
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Eq  (66)  is  the  fractional  order  state  equations  with  optisal  feedback 
control.  Application  of  this  control  law  will  sake  possible  the 
controlling  of  fractional  derivatives  of  the  structural  response. 

3.2  Fractional  Order  State  Feedback  Control  of  Viscoelastically  Daaped 
Structures 


The  next  step  in  the  study  is  to  forsuiate  the  fractional  order 
state  feedback  control  law  for  the  control  of  structures  that  are 
viscoelastically  daaped.  The  development  process  follows  closely  to  the 
procedure  used  in  the  development  of  the  fractional  order  state  feedback 
control  law  for  the  control  of  undamped  structures.  Minor  modification 
was  made  in  the  formulation  of  the  finite  element  model  of  the 
viscoelastically  damped  structures  which  included  the  fractional 
derivative  modeling  of  the  viscoelastic  damping  material. 

Consider  the  three  parameter  fractional  derivative  model  for  a 
viscoelastic  material  undergoing  uniform,  uniaxial  deformation 

o(t)  =  (Eq  +  BtD^)e(t)  (67) 

where  E0,Elt  and  0  are  the  material  parameters.  Applying  this  model  to 
the  formulation  of  the  finite  element  equations  of  notion  for  a 
viscoelastically  daaped  structure  under  control  and  the  equations  take 
the  fora 

lMl{w(t))  +  [KvHl^w(t))  ♦  [Kg]{w(t))  »  tD](u(t))  (68) 
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where 


[Ky]  =  visco-st iff ness  matrix 
[Ke1  =  elastic-stiffness  aatrix 

In  general  the  visco-stiffness  aatrix  is  not  a  linear  combination  of  the 
■ass  and  elastic-stiffness  matrices,  and  the  equations  of  motion,  Eq 
(68),  can  not  be  transformed  into  the  traditional  first  order  state 
equations,  Eq  (1).  In  order  to  apply  optimal  feedback  control  to  this 
viscoelastically  damped  structure,  it  is  necessary  to  develop  an 
equivalent  damping  matrix  and  an  equivalent  stiffness  matrix  for  the 
system,  thus  the  equations  of  notion  may  then  be  posed  in  the  familiar 
form 

CM]{w(t))  +  [C*)Mt)>  +  [K*]{w(t) )  *  [DHu(t))  (69) 

and  Eq  (69)  could  lead  to  the  formulation  of  the  necessary  state 
equations  for  structural  control. 

To  develop  the  equivalent  damping  and  stiffness  matrices,  one  must 
first  solve  for  the  resonance  response  of  the  viscoelastically  damped 
structure.  Posing  the  equations  of  notion  in  the  fractional  format 


Dy[M] 


-  [M] 


♦  D2'22  ♦  D2'3*  ♦  ■—  ♦  D°](.(t)| 

[d2-»  .  D2-2»  .  D2*3*  .  -♦»»]{«<»)  I 

,  D»[Ky]  .  DP'31  .  -  .  D°] 

(w(t)) 

-  [Kvlp-»  *  D?-2»  .  D?-3>  *  -  , 

Dr]{w(t)} 

+  [KE]D°{w(t)>  =  {0} 


(70) 


where  7  =  1/a  and  a  is  equal  to  the  denoainator  of  0.  The  aodal 
displacements,  modal  velocities,  and  aodal  accelerations  of  the 
structure  can  then  be  deterained  by  solving  Eq  (70)  using  the  solution 
aethod  of  the  expanded  equations  of  notion  presented  in  Eq  (32)-Eq  (36). 
The  expanded  equations  of  notion  for  Eq  (70)  are  written  as 


DT[M]{w(t)}  +  [K]{w(t)l  *  (0) 


(71) 


where 
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(72) 
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.  “[Kyi  ["1 

[01 
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and  (w(t)}  is  the  state  vector  defined  in  Eq  (35).  Once  the  expanded 
eigenvalues  and  eigenvectors  are  determined  by  solving  Eq  (71),  one  say 
proceed  to  establish  the  modal  matrix,  [♦] ,  for  the  uncontrolled  system 
as  described  in  Eq  (39)-Eq  (41)  and  the  nodal  displacements,  velocities, 
and  accelerations  can  be  extracted  from  the  nodal  matrix  and  the 
expanded  eigenvalues  on  the  principal  sheet  of  the  Rienann  surface. 

Given  that  the  equivalent  damping  and  stiffness  matrices  must 
generate  the  same  structural  response  as  the  original  system,  one  may 
solve  for  the  matrices  provided  that  the  resonant  structural  response  of 
the  original  system  is  known.  The  homogeneous  equations  of  motion  of  the 
equivalent  system  may  be  written  as 

[M]D2{w(t)}  +  [C*]D{w( t) }  ♦  [K*l{w(t)>  *  {0}  (74) 

Applying  the  expansion  theorem  for  structural  modal  analysis  (24:283)  to 
Eq  (74)  and 


[M][*j{n)  ♦  [c*][*vHn>  ♦  [K*][OdHn>  *  (o> 


(75) 


Further  simplification  produces 
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(76) 


(M][*J  +  [C*J[*v]  +  [K*][*d]  =  {0} 

where  [0^],  [$y],  and  [$dl  are  the  nodal  acceleration,  velocity,  and 
displaceaent  aatrices  respectively.  These  aatrices  consist  of  the 
eleaents  of  the  aodal  aatrix  associated  with  the  known  structural  nodal 
acceleration,  velocity,  and  displaceaent.  Since  the  eleaents  of  the 
aodal  aatrix  are  coaplex,  Eq  (76)  Bay  be  redefined  as  two  equations 

[M]Re[*  ]  +  [C*]Re[4>  ]  +  [K*]Re[*J  *  {0}  (77) 

•  v  a 

[M]In[*J  +  [C*]IbI$v]  ♦  [K*]Ia[*d]  =  {0}  (78) 

where  Re[^J,  Re[<M ,  and  Re[^]  represent  the  real  parts  and  InOJ, 
Ia[$  ],  and  Ia[$.]  represent  the  iaaginary  parts  of  [♦  J,  [♦  J  and  [♦ ] 
respectively.  Eq  (77)  and  Eq  (78)  Bay  now  be  solved  for  [C*]  and  [K*] 
where 


[C*]  =  -[Mj^ReCfJ  -  Ia[^)  lB[4d]‘t  Re[^d]j 

•[Be[#vJ-Ia[#T]-lBC#d]'l*Re[#d])'t  (79) 

[K*]  *  -([MWa[*J  +  [C*]  lB[^]]lB[^dJ_1  (80) 

After  the  equivalent  daaping  and  stiffness  aatrices  have  been 
deterained  for  the  systea,  the  state  equations  for  the  equivalent  systea, 
Eq  (69),  can  be  developed  as 
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X  *  Ax  +  Bli 


(81) 


where 


[0]  [I] 

.-wrV]  "[m]"1[c*j  . 


(82) 


[0] 

.  (M1_1[D1 


(83) 


The  optimal  gain  aatrix,  G,  for  feedback  control  can  be  deterained  in 
the  saae  Banner  aa  described  in  Section  2.1.  Incorporating  the  gain 
aatrix  into  the  systea  aodel,  and  the  finite  eleaent  aodel  of  the 
controlled  viscoelastically  daaped  structure  takes  the  fora 


[M]{w(t)>  4  ([C*]  +  [D][G2]){w(t)} 

♦  ([K*]  +  [D][G1]){w(t)}  =  (0)  (84) 

To  deteraine  the  state  transforaation  aatrix  for  the  systea,  one 
should  first  pose  Eq  (84)  in  the  fractional  foraat 

AmiJd2-*  ♦  D2-**  4  D2-3"  +  -  4  D°]{w(t)> 

-  [M][D2-a  4  D2-2a  4  D2"30  4  --  4  Da]{w(t)) 

4  Da([C*]  4  [D][G2] )  [d1-<*  4  D'-2*  4  D1'3*  4  4  D° j  { w ( t )  } 

-  (tc*l  4  [DltGg] ) 4  Dl'2a  4  D1_3a  4  •—  4  Daj{w(t)} 

4  ( [K*l  4  [Dj[Gtl)D°(w(t)}  »  (0)  (57) 
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The  nodal  transformation  aatrix  aay  be  formulated  by  solving  the 
fractional  order  equations  of  notion  in  Eq  (85),  using  the  solution 
nethod  of  expanded  equations  of  notion  in  the  manner  described  in  Eq 
(32)-Eq  (36).  The  expanded  equations  of  notion  for  Eq  (85)  are 


'V  -v 


Da[M]{w(t) >  +  [Kl{w(t)>  =  (0) 


(86) 


where 


[M]  = 


[K]  = 
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[0] 
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tl+[Dl[G2] 

[01 
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[C*MDHQJ 

[01 

[01 

[M] 

[01 

[c*]+[Dl(G2l 

[01 

[01 

[0] 

[0J 

— 

[01 

-[M] 

[01 

[01 

[01 

— 

-[Ml 

. 

[0] 

j 

| 

!  1  ! 

j  - 

[c*i-[d][g2i 

j 

[01 

-[Ml 

- . [C*J-[D][G2J  [01 

[0] 

[0] 

-[M] 

[01 

—  -[C*]-[D][G21 

[01 

[0] 

[01 

[01 

[0] 

—  .  [01 

[01 

[K*MD][G 

(87) 


(88) 


Following  the  procedure  identified  in  Eq  (39)-Eq  (48),  the  state 
transfornation  natrix,  ♦,  can  be  deternined  fron  the  systea  nodal 
natrix,  and  the  control  law  for  the  feedback  of  fractional  order  states 
nay  be  written  as 


33 


x  =  Ax  -  BG$x 


(89) 


The  control  law  of  Eq  (89)  for  the  control  of  viscoelastically  daaped 
structures  remains  unchanged  from  the  control  law  for  the  control  of 
undaaped  structures,  only  the  values  of  the  optimal  gain  matrix  and  the 
state  transformation  are  different. 
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IV.  Numerical  Approxiaation  of  D^[ezt] 


The  purpose  of  this  section  is  to  develop  nuaerical  approximation 
of  the  fractional  derivative  of  the  complex  exponential  function 

F(t)  =  D®[ezt]  (90) 

for  0  <  <x  <  1  and  s  is  a  complex  number.  The  nuaerical  approxiaation  of 
Eq  (90)  provides  the  foundation  for  the  investigation  of  general 
behavior  of  structures  that  are  viscoelastically  damped  and  structures 
that  are  controlled  using  fractional  order  state  feedback  control 
theory.  It  can  be  shown  herein  that  the  behavior  of  the  system  is 
strongly  influenced  by  the  poles  on  the  principal  sheet  of  the  Riemann 
surface  and  the  residues  of  these  poles  provide  good  approximation  for 
the  general  system  response. 

In  this  section  two  nuaerical  approxiaation  expressions  for  the 
fractional  derivative  of  the  complex  exponential  function  are  developed. 
These  expressions  are  used  to  approximate  the  short  range  and 
intermediate  range  behavior  of  the  function  in  Eq  (90).  The  expressions 
for  the  errors  associated  with  these  approximations  are  also  developed 
in  the  process.  Finally  these  approximations  are  used  to  generate 
estimations  for  the  fractional  derivative  of  the  generalized  sinusoidal 
function 


G(t)  *  D®[eli  +  ezt] 


(91) 
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where  ezt  is  the  cosplex  conjugate  of  ezt.  Recalling  fros  vibration 
analysis  (24:10)  that  the  general  response  of  an  underdasped  system  is 


X(t)  =  A^e*1  ♦  A2ezt 


(92) 


where 


z  =  -C«n  +  i«d  =  (-C  +  i(l  -  CV'2  )  s  3  wn 
«  =  -C«n  -  i<*d  =  (-C  -  id  -  C2)1'2  )  «n  -  i  «n 


(93a) 

(93b) 


and 

3  *  (-C  ♦  id  -  CV'2  )  (94a) 

I  ■  (-C  -  id  "  CV'2  )  (94b) 

where  C  is  the  damping  ratio  of  the  systes.  and  wd  are  the  natural 
frequency  of  the  undamped  system  and  the  resonant  frequency  of  the 
damped  system  respectively.  Thus  Eq  (91)  can  provide  a  general 
description  of  the  fractional  derivatives  of  the  response  of  an 
underdamped  structural  mode. 

Before  proceeding  with  the  development  of  the  approximation  for  the 
fractional  derivative  of  the  complex  exponential  function,  it  is  useful 
to  pose  the  exponential  function  in  terms  of  the  Maclaurin  series 
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The  last  step  in  calculating  the  derivative  in  Eq  (90)  is  to  calculate  the 
inverse  Laplace  transfora  of  Eq  (98) 
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(99) 


(100) 


Let  a  =  n  -  a  in  Eq  (100)  and  apply  the  inverse  transfora  to  Eq  (99) 
then 


(n-a) 


_ar  zt,  \  n  t 

D.te  1  *  1  *  r<n--«  ♦  1) 


n»0 


(101) 


Eq  (101)  is  the  fractional  derivative  of  a  complex  exponential  function. 
Express  Eq  (101)  in  terns  of  its  Nth  partial  sun  and  its  remainder 


_ar  *t. 
Dje  ] 


■2* 

n»0 


(n-a) 


T(n  -  a  ♦  1) 


,  ?  ■  t(n-g) 

2  *  l"(n  -  a  +  1) 


(102) 


where  the  first  part  of  Eq  (102)  is  the  Nth  partial  sum  of  the  series 
and  the  second  part  of  the  equation  is  the  remainder  of  the  series  after 
tern  N. 

Redefine  Eq  (102)  in  terms  of  dimensionless  tine,  0(t),  where 


0(t)  *  «nt 


(103) 
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Eq  (102)  aay  then  be  written  as 


D«te3«(t), 


K 


><n-0O 


a> 


e 

fyn  -  a 


+  1) 


n«0 


K 


n-H  +  1 


0 


(n-a) 


r(n  -  a  ♦  1) 


(104) 


The  general  approxiaation  of  the  fractional  derivative  of  the  coaplex 
exponential  function  can  be  expressed  as  the  Nth  partial  aua  of  the 
series 


«  *  n  . (n-a) 

_ar  zt,  ^  V  n  t 

®t[e  1  *  1  ‘  f(n  -  a  ♦  1 

n»0 


) 


*  yu«  3"_e(.nIa.— 

Z  n  "  T(n  -  a  +  1 

n»0 


) 


(105) 


and  the  reaainder  of  the  series  is  the  error,  EN,  of  the  approxiaation 


OD 


g  .  X  _  t’-"1 

E«  Z  *  f(n  -  «  ♦ 


1) 


(n-a) 


I  Wn  3  T(n  -  a  ♦  1) 


n«M*l 


(106) 


Notice  that  the  error  function  EN  increases  without  bound  as  0(t) 
increases,  and  the  nuaber  of  teras  required  to  aaintain  a  prescribed 
degree  of  accuracy  in  the  solution  increases  rapidly  as  0(t)  increases. 
Therefore  Eq  (105)  is  only  valid  for  the  short  range  approxiaation  of 
the  fractional  derivative  of  the  coaplex  exponential  function. 
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To  better  understand  the  general  behavior  of  the  function,  Eq  (90), 
it  is  necessary  to  develop  an  expression  that  can  better  portray  the 
overall  behavior  of  the  function.  The  process  begins  by  taking  the 
Laplace  transform  of  Eq  (90) 


e[D®[eZt]]  *  sa-2[ezt] 


(107) 


The  Laplace  transform  of  the  exponential  function  is  determined  froa 


or  zt,  f  zt  - •  t  1 

P[e  ]  =  I  e  e  dt  *  - 

J  8*8 


(108) 


Substitute  Eq  (108)  into  Eq  (107)  and  the  result  is 


8  -  Z 


(109) 


To  determine  the  fractional  derivative  of  the  complex  exponential 
function,  one  needs  to  calculate  the  inverse  Laplace  transform  of  Eq 
(109).  Following  the  definition  of  inverse  Laplace  transform  integral 


e’t[F(s)]  *  J  [F(s)J  e-t  ds 
y-iw 


(110) 


the  inverse  transform  of  Eq  (109)  may  be  expressed  as 


40 


y+i  u 


*•“> .  J  (rh-)  e“  * 


Y-  iu 


(111) 


and  Eq  (111)  con  be  evaluated  using  the  residue  theorem  from  the 
calculus  of  coaplex  variables.  The  closed  contour  of  integration  for 
the  inverse  transform,  in  conjunction  with  the  residue  theorem ,  is  shown 
in  Figure  1. 

The  integration  contour  has  six  segments  and  the  direction  of 

I 

integration  is  indicated  by  the  arrows.  Segments  3,  4,  and  5  are 
included  in  the  contour  because  the  branch  cut  of  the  function  s1/a  is 
token  along  the  negative  real  axis  and  the  branch  point  of  the  sane 
function  is  taken  at  the  origin  of  the  s  plane.  The  inverse  transfora, 
Eq  (111),  is  represented  by  segment  1  of  the  contour. 
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The  residue  theorem  states  that  the  integral  along  the  contour, 
divided  by  2Jti,  i»  equal  to  the  sue  of  the  residues  of  the  poles  of  the 
integrand  enclosed  by  the  contour  (22:865).  By  including  the  inverse 
transfora  in  the  integration  contour  and  applying  the  residue  theorea 
froa  coaplex  calculus,  the  inverse  transfora,  Eq  (111),  Bay  be  written 
as 


y+iw 


zfe  I  (r^i)  e“  d< *  "~2kI  I  I  (r?i) e’1  d* + 1  bj  (ll2> 


7-iw 


k-2  k 


where  the  first  suaaation  represents  the  sua  of  the  integrals  of  the 
other  five  segaents  and  the  second  suaaation  represents  the  sua  of  the 
residues.  The  residues,  ,  are  evaluated  by  conventional  techniques 

b.  =  lia  (s  -  IJlPtsHe'1  (113) 


provided  the  liait  exists.  Under  the  current  analysis,  there  is  one 
pole  at  Xj  *  z  and  the  residue  of  this  pole  enclosed  by  the  contour  is 

bj  =  lia  (s-s)  f-  *  |  e#t  *  zae2t  (114) 

It  has  been  shown  (5:95-99)  that  the  integrals  along  segaents  2,  4, 
and  6  are  equal  to  zero  when  the  radii  of  2  and  6  are  increased 
indefinitely  and  the  radius  of  4  is  decreased  indefinitely.  The 
integrals  along  segaents  3  and  5  are  not  equal  to  zero  and  the  sua  of 
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the  integrals  are  not  equal  to  zero  either.  The  expression  for  the 
inverse  transform  given  by  Eq  (112)  is  now 


ir+io 


1 

2ni 


•  t 


ds  + 


«  zt 
z  e 


(115) 


T-iW 


k-3,5  b 


and  the  inverse  transfors  is  now  the  sus  of  the  integrals  along  the 
branch  cut,  segments  3  and  5,  plus  the  residue.  The  sus  of  the 
integrals  can  be  deterained  as  follow 


-isr  l  J  (rh)  •** d< 


k-3,5  k 
/  O 
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2ni 


(re1*)*  e(reiKt)jn 
i;  (reiK-z) 


CD 
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in. a 
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e-  dr  ♦  f-lg  „<re  *>, 

{  (re"*-!) 
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~JL  f  - 

2ni  [ 


CD 


ina  ~ma 
e  +  e 


)Jff 


a  -rt 


dr 
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<9 


*  -sir  ( Jffrdr 


a  -rt 


=  -T*in(°*>  !tTT  dr 


(116) 


Applying  the  reflection  formula  of  gamma  function  (27:256)*  where 


r<*)ra  -  >)  ■  -*r(-s)r(s)  * 


(117) 


and  Eq  (116)  can  be  written  as 


i 
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) 


1 

2ni 


2  ffc£) 


•  t 

e 


ds 


tVirll  -  ~)  J-z(f:'rT») dr  (118) 


Recalling  the  geometric  series  where 


n*0 


1  +  X  +  X 


(119) 


Let  x  =  -r/z  in  Eq  (119)  and  apply  Watson’s  lenna  (30:21-24)  to  Eq  (118) 


1ST  1  I  (rh)  «‘>d>  -  r(g)f(i  - »)  H  ■  <12°> 


b-3,5  k 


nmO 


Now  substitute  Eq  (120)  into  Eq  (115)  and  the  fractional  derivative  of 
the  coapiex  exponential  function  can  be  expressed  as 


CD 


„ar  It,  ..  o  *t  .  1  VI  r(a+n+l)  .  .  _ 

“J*  >""e  ♦  r<«)fu  - «)  ZT  (121) 


Eq  (121)  is  the  asynptotic  expansion  representation  of  the  saae  function 
in  Eq  (90).  Redefine  Eq  (121)  in  teras  of  diaensionless  tine,  0(t), 
defined  in  Eq  (103),  and  Eq  (121)  can  be  written  as 


CD 


a*!> 


°>38,t,>  -  •:  >°  ♦  Twa  -  «>  i  *  -  -  ««*, 


n»0 


I 
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The  general  behavior  of  the  fractional  derivative  nay  be 


approx i sated  aa 


_<*,  z  t .  a  zt 
D  [e  ]  *  z  e 


(123) 


for  the  intersediate  range  of  tine  and  the  sunaation  ters  represents  the 
reaainder,  R  ,  of  the  approxiaation 
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1 

CD 

_  \± 

r(a+n-H) 

'  nooni 

-  «)  Z  z 

n«0 

o+n*i 

i 

CD 

_  V  _ 

n  r(a+n+l) 

-  r(o)T( i 

-or )  L 

n*0 

8  8“*"*' 

(124) 


The  reaainder  function  R  now  decreases  as  0(t)  increases  for  this 

n 

function.  Because  the  exponential  tern  decays  faster  than  the  suaaation 
tern,  the  long  range  behavior  of  the  function  is  influenced  more  by  the 
reaainder  tern.  The  effect  of  the  reaainder  is  even  aore  pronounced  at 
low  frequency  and  in  short  tine  range.  Thus  Eq  (123)  provides  only  a 
representation  of  the  interaediate  range  behavior  of  the  function  under 
investigation. 

Given  Eq  (105)  and  Eq  (123),  one  aay  proceed  with  the  evaluation  of 
the  approxiaation  of  fractional  derivative  of  the  generalized  sinusoidal 
function  defined  by  Eq  (91).  Since  z  and  z  are  coaplex  conjugates,  the 
short  range  approxiaation  of  Eq  (91)  for  systeas  with  zero  initial 
displaceaent  can  be  written  as 
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(125) 


D®[ezt  +  ezt]  *  2-^  Ia(z") 

n»0 


(n-a) 


f(n  -  a  ;  1) 


and  the  approxiaation  in  Eq  (125)  expressed  in  diaensionless  tine  0(t) 
and  normalized  by  2  is 

TV 


,  -a  _ar  zt 
"n 


+ 


a  Q(n-a) 

3  p  _ 

f(n  -  a  +  1) 


n  «  0 


) 


(126) 


In  the  saae  manner,  the  intermediate  range  approximation  of  Eq  (91)  can 
be  posed  as 


_ar  zt  ,  zt. 

Dt[e  ♦  e  ] 


*  24Im(zaeZt) 


(127) 


and  the  approxiaation  of  Eq  (127)  expressed  in  dimensionless  time  0( t) 


a 


and  normalized  by  2-u  is 


1  nar  zt  i  ztl  -  r  i*a  B®, 

—  <*>n  Dje  «•  e  ]  "  Im(B  e  ) 


(128) 


The  normalized  approxiaation  of  the  short  range  and  intermediate  range 
behavior  for  the  fractional  derivative  of  the  response  of  an  underdaaped 
system,  Eq  (126)  and  Eq  (128)  respectively,  with  a  =  0.5,  C  s  0.1, 

*  x  and  N  *  60  are  illustrated  in  Figure  2.  Notice  that  the 
intermediate  range  approxiaation  converges  onto  the  exact  behavior 
expression  in  less  than  one  cycle. 
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Figure  2.  Approximation  of  Fractional  Derivative  of 
Generalized  Sinusoidal  Function 


Based  on  the  analysis  just  presented,  one  say  conclude  that  the 
systea  response  represented  by  the  intermediate  range  function  is 
a  relatively  accurate  representation  of  the  actual  response  of  the 
system  studied,  and  the  response  is  strongly  influenced  by  the 
residues  of  the  poles  of  the  systea.  This  relationship  forms  the 
basis  of  the  fractional  order  state  feedback  control  theory  formulation 
in  Chapter  III,  where  only  the  poles  on  the  principal  sheet  of  the 
Rieaann  surface  are  retained  for  the  development  and  formulation  of  the 
state  transformation  matrix.  Following  this  analysis,  one  nay  formulate 
the  fractional  order  state  feedback  control  force  input  based  on  the 
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system  response  and  the  eigenvalues  and  eigenvectors  associated  with  the 
poles  on  the  principal  sheet  of  the  Rienann  surface  alone,  and  still 
Maintain  equivalent  control  authority  as  the  traditional  state  feedback 
control  with  fine  control  accuracy. 
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V.  Example  Probleas 

To  deaonstrate  the  application  of  the  fractional  order  feedback 
control  theory,  the  finite  eleaent  analysis  of  two  siaple,  controlled, 
structures  are  presented  in  this  section.  The  nuaerical  solution  of  the 
exaaples  are  coaputed  using  aatrix  aanipulation  routines  in  the  prograa 
PRO-MATLAB  (26).  The  first  exaaple  is  a  siaple,  hypothetical,  two 
degree-of-freedon  spring-aass  systea  undergoing  longitudinal  vibration. 
Both  active  closed-loop  feedback  control  and  passive  viscoelastic 
daaping  coupled  with  active  feedback  control  are  exaained  in  this 
exaaple.  The  second  exaaple  is  part  of  a  structural  control  experiaent 
currently  being  developed  at  AF1T  to  deaonstrate  the  applicability  of 
the  fractional  order  state  feedback  control  theory  in  controlling  large 
flexible  space  structures.  The  objective  of  this  exaaple  is  to  develop 
the  optiaal  gain  and  state  transforaation  aatrices  for  the  control  of 
the  aodified  Advanced  Beaa  Experiaent  (ABE)  (25)  using  fractional  order 
state  feedback. 

Exaaple  Problea  1 

In  this  first  exaaple,  the  application  of  fractional  order  state 
feedback  control  theory  in  structural  control  will  be  deaonstrated, 
using  nuaerical  analysis  of  a  hypothetical  structure.  The  first  part  of 
the  exaaple  will  illustrate  the  foraulation  of  the  expanded  equations  of 
notion,  the  deteraination  of  the  resonant  frequencies  and  node  shapes, 
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and  the  coaposition  of  the  state  transformation  matrix  for  an  open-loop 
system. 


Consider  a  two  degree-of-freedoa  system,  shown  in  Figure  3,  in 
longitudinal  free  vibration.  The  mass  and  stiffness  matrices  for  the 
system  are  defined  as 


[MJ  = 


[K]  * 


0.5  0 

0  0.5 

6  -2 
-2  6 


The  general  form  of  the  equations  of  notion  for  the  homogeneous  system 
is 

(Ml{w(t)}  +  [K]{w(t)}  *  {0}  (129) 

where  (w(t)}  is  the  displacement  vector  defined  as 
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{w(tn 


(130) 


The  equations  of  aotion  for  the  system  under  free  vibration  aay  be  posed 
as 


'  0.5 

0 

6 

-2  ' 

< 

+ 

< 

.  0 

0.5  . 

N 

.  -2 

6  . 

N 

and  the  eigenvalues  and  eigenvectors  of  the  open  system  can  be 
deterained  by  the  EI6  coaaand  of  PRO-MATLAB  (26:3.35).  These  values 
correspond  to  the  first  two  aodes  of  vibration.  The  natural  frequencies 
of  the  8ystea  are 

ut  =  2.8284 

«  *  4.0000 

£ 

and  the  node  shapes  of  the  system  are 

Mz  *  {-! } 

The  general  problea  can  also  be  solved  by  first  posing  the 
equations  of  notion  in  a  fractional  order  fornat  and  solving  the 
fractional  order  equations  using  the  solution  technique  of  the  expanded 
equations  of  notion.  Assuming  a  s  1/2  and  applying  Eq  (25)  in  the 
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equation  formulation,  Eq  (129)  aay  be  written  as 


D1/2CM][d3'2  +  D1  Dt/2  ♦  D°J{w(t)} 

-  CM]|d3'2  +  D1  +  D1/2]{w(t)> 

+  [K]D°{w(t) )  =  (F(t) }  (132) 


To  solve  for  the  eigenvalues  and  eigenvectors  of  Eq  (132),  follow  the 
solution  aethod  of  expanded  equations  of  notion  described  in  Eq  (32)-Eq 
(36)  and  arrange  Eq  (132)  in  the  expanded  forsat 


D,/2[M](w(t)}  +  (K](w(t) )  ={F(t>  > 


where 


[0] 

[0] 

[0] 

CM] 

CM]  = 

[0] 

[0] 

CM] 

CO] 

[01 

CM] 

[0] 

CO] 

CM] 

[0] 

[0] 

co] 

’  [0] 

[0] 

-CM] 

[0] 

IKJ  » 

[0] 

-CM] 

CO] 

co] 

-CM] 

[0] 

[0] 

co] 

[0] 

[0] 

[0] 

CK] 

(133) 


(134) 


(135) 
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{w( t ) >  = 


> 


(136) 


'  D3/2{w(t)) 

D1  <w(t) ) 
D1/Z{w(t)) 

„  (w(t)}  J 


The  eigenvalues  and  eigenvectors  of  the  expanded  equations  can  then  be 
solved  using  the  EIG  consand  of  PRO-MATLAB  and  the  eigenvalues  are 

X  a  -  1.4142  ±  1.41421 

X  a  -  1.1892  ±  1.18921 

X  a  1.4142  ±  1.41421 

X„  »  1.1892  ±  1.18921 

7.9 

The  eigenvalues  with  positive  real  parts  say  be  sapped  onto  the 

1/2 

principal  sheet  of  the  Rieaann  surface  for  the  function.  X5  =  s 
consisting  of  two  sheets  (23:902).  These  eigenvalues  and  the 
corresponding  eigenvectors  describe  the  sinusoidal  notion  of  the 
structure  and  fors  the  basis  for  the  developsent  of  the  state 
transformation  matrix  as  described  in  Chapter  III  of  this  thesis.  The 
natural  frequencies  of  the  system  can  be  determined  by  squaring  the 
fifth,  sixth,  seventh,  and  eighth  eigenvalues.  This  process  produces 
two  conjugate  pairs 

X?  *  0.0000  ±  4.00001 

5,0 

X*  »  0.0000  ±  2.82841 

7,8 
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The  imaginary  parts  of  each  conjugate  pair  correspond  to  a  natural 
frequency  of  the  systea.  The  zero  in  the  real  part  indicates  that  no 
daaping  is  present  in  the  structure. 


The  eigenvectors  associated  with  the  eigenvalues  identifying 


systea  resonant  behavior  (fifth,  sixth,  seventh,  and  eighth)  also  occur 


in  coaplex  conjugate  pairs.  They  nay  be  posed  as  a  nodal  natrix  of  the 
fora,  Eq  (39),  and  the  aodal  aatrix  for  the  systea  is 


[<J>]  = 


1.0000 

1.0000 

1.0000 

1.0000 

-1.0000 

-1.0000 

1.0000 

1.0000 

0.3536 

- 

0.35361 

0.3536 

♦ 

0.35361 

0  4204  -  0.4204i 

0.4204 

+ 

0. 4204i 

-0.3536 

+ 

0.35361 

-0.3536 

- 

0.35361 

0.4202  -  0.42021 

0.4204 

+ 

0.42041 

0.0000 

- 

0.25001 

0.0000 

+ 

0.2500i 

0.0000  -  0. 3536i 

0.0000 

+ 

0.35361 

0.0000 

+ 

0.25001 

0.0000 

- 

0.25001 

0.0000  -  0.35361 

0.0000 

+ 

0.3536i 

-0.0884 

- 

0.08841 

-0.0884 

♦ 

0.08841 

-0.1487  -  0. 1487i 

-0.1487 

♦ 

0. 1487i 

0.0884 

+ 

0.08841 

0.0884 

- 

0.0884i 

-0.1487  -  0.14871 

-0.1487 

+ 

0.14871 

The  seventh  and  eighth  rows  of  the  aatrix  represent  the  node  shape  of 

the  systea.  The  fifth  and  sixth  rows  are  associated  with  the  node  of 

t  /  2  t 

the  D  state.  The  third  and  fourth  rows  represent  the  node  of  the  D 

»  3/2 

state  and  the  first  and  second  rows  represent  the  node  of  the  D  state 

of  the  systea. 


For  this  open  systea,  the  state  equations  can  be  written  as 


k  *  Ag 


(137) 
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where 


X  = 


A 


> 


(138) 


Recalling  froa  the  expansion  theorem  in  structural  aodal  analysis  that 
the  state  vector  nay  be  posed  as  the  product  of  an  associated  nodal 
natrix  and  noraal  coordinates,  the  state  vector  in  Eq  (138)  can  be 
written  as 


X  *  [*A]{n> 


where 


-0.0884  -  0.08841 
0.0884  -i-  0.08841 
0.3536  -  0.35361 
-0.3536  +  0. 3536i 


-0.0884  +  0.08841 
0.0884  -  0.0884i 
0.3536  *  0. 3536i 
-0.3536  -  0. 3536i 


-0.1487  -  0.14871 
-0.1487  -  0. 1487 i 
0  4204  -  0.42041 
0.4202  -  0.4202i 


(42) 


-0.1487  ♦  0. 1487i 
-0.1487  +  0. 1487i 
0.4204  *  0. 4204i 
0.4204  +  0. 4204i 


for  the  present  systen.  The  row  eleaents  of  [#A]  are  those  row  eleaents 
of  [♦]  that  correspond  to  the  states  of  the  state  vector  &.  In  a 
similar  Banner,  the  fractional  state  vector  nay  be  written  as 


i  *  [tg]{nf 


(44) 


i 
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For  the  application  of  one  fractional  state  observer  placed  at  y  can 
be  defined  as 


1  = 


V'2 


(139) 


and 


1.0000 

0.3536  -  0.35361 
0.0000  -  0. 2500i 
-0.0884  -  0.0884i 


1.0000 

0.3536  +  0. 3536i 
0.0000  +  0.25001 
-0.0884  +  0.08841 


1.0000 

0  4204  -  0.4204i 
0.0000  -  0.35361 
-0.1487  -  0. 1487i 


1.0000 

0.4204  +  0. 4204i 
0.0000  ♦  0. 3536i 
-0.1487  +  0. 1487i 


for  this  systes.  Here  the  row  elesents  of  [♦g]  are  those  row  elesents 
of  [0]  associated  with  the  state  elesents  identified  in  y. 

To  deteraine  the  state  transformation  matrix,  <P,  between  the 
traditional  state  vector  and  the  fractional  state  vector*  apply  Eq  (48) 
and 


<p  3 


1.3017 

0 

4.3784 


0  0  1 

-3.3889  4.3784  0.1716 

10  0 

-22.9693  44.6249  -38.3410 
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The  state  transformation  matrix  allows  one  to  determine  the  overall 
structural  response  by  observing  only  the  fractional  and  integer  states 
at  m^  and  this  relationship  can  be  illustrated  as 

X  ■  (49) 

Eq  (49)  can  be  defined  for  this  example  as 
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In  the  second  part  of  this  exaaple,  a  control  force  is  applied  at 
■t  for  the  purpose  of  implementing  feedback  control  to  the  system.  The 
system  set  up  is  shown  in  Figure  4. 


u. 


2K 

ww 

1/4 


WW 

L/2 


2K 

W/W 

1/4 


Figure  4.  Feedback  Controlled  2  OOF  System 


The  equations  of  notion  for  the  system  are 


[M]{w(t) }  +  [K] { w( t ) }  =  [D]{u(t) } 


(141) 


and 


[D]  = 


(142) 


The  one  entry  on  the  diagonal  indicates  that  the  location  of  control 
force  is  at  nt.  The  state  model  of  the  equations  of  notion  may  be  put  in 
the  form 


X  *  A*  ♦  Bu 


(1) 
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and  the  elements  of  the  state  aodel  are  defined  in  Eq  (51),  (52),  (53), 
and  (54). 

The  optiaal  control  for  the  system  is  deterained  by  ainiaizing  the 
quadratic  performance  index 


J  =  yJ[x(t)TQx(t)  +  uT(t)Ru(t)]dt  (143) 

o 


where 


Q  * 

[K] 

[0] 

.[01 

[Ml. 

D  • 

100 

0 

K  — 

0 

100 

(144) 


(145) 


Notice  that  the  first  term  of  the  integral  represents  the  total  energy 
of  the  system  and  the  second  integral  tern  represents  the  control 
energy.  The  optiaal  control  for  the  system  is 


a  *  -  Gx 


(7) 


where  G  is  the  optimal  control  gain  aatrix  and  it  is  defined  as 

G  *  R_1BS  (4) 

and  3  is  the  solution  of  the  aatrix  Biccati  equation.  For  the  current 
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systea  the  optiaal  gain  aatrix  is 


G 


[G, 


Gz] 


0.0149  -0.0298 

0.1412 

0.0002 

0  0 

0 

0 

Only  the  first  row  of  the  gain  aatrix  is  populated  because  only  one 
control  force  is  applied  to  the  systea. 

Incorporating  the  optiaal  feedback  control  into  the  systea,  the 
equations  of  notion  nay  now  be  written  in  a  homogeneous  fora 

(M]{w(t) }  ♦  [D][G2]{w(t)}  +  ([K]  +  [DHGt]){w(t)}  =  {0}  (56) 

The  systea  response  can  be  determined  by  first  posing  Eq  (56)  in  the 
fractional  format  and  then  solving  the  problem  using  the  solution  method 
of  the  expanded  equations  of  notion  described  in  Eq  (32)-Eq  (36). 

For  a  =  1/2,  the  fractional  order  equations  of  notion  are 

D^tMl^  +  D1  ♦D,/2  +  D°]{w(t» 

-  [M]|d3/2  +  D1  +  D1/ZJ{w(t)} 

+  D1'2[D][G2][d1'Z  +  D°]{w(t>} 

-  [Dl[G2]D1/2{w(t)} 

+  UK]  ♦  (D][G1])D°{w(t)}  *  {0}  (146) 

and  the  expanded  equations  of  notion  are 
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f  /  2  *V  «<V  -V  -V 

D1/^[M]{w(t)}  +  (K]{w(t) }  =  {0) 


where 


[MJ  = 


[K]  = 


[0] 

to] 

to] 

tM] 

[0] 

to] 

tM] 

to] 

[01 

tMl 

10] 

tD]  [G 

[M] 

to] 

tD] tG2]  t0] 

[0] 

[0] 

-tM] 

[0] 

~[M] 

[0] 

-[M] 

[0] 

- 

tD]tG2] 

[0] 

to] 

[01 

{w(t} ]  = 


D3'2{w(t)} 
D1  <w(t)} 
D1/2{w(t)} 


There  are  four  pairs  of  cosplex  conjugate  eigenvalues  to  Eq  (147) 
and  they  are 


S  -  1.4018  ±  1 . 426 7 i 


a  -  1.1745  ±  1 . 204 1 i 
=  1.4018  ±  1 . 4267 i 


1.1745  ±  1 . 204 1 i 


The  eigenvalues  X  and  X„  correspond  to  the  node  two  and  node  one 

0,0  7,0 

behaviors  of  the  systea.  Squaring  the  eigenvalues  generates  the 
respective  frequencies  for  the  nodes 


X"  =  -  0.0706  ±  4.0000i 
o ,  o 

X*  a  =  -  0.0706  ±  2.8284i 

7 , 0 


(150) 


where  the  iaaginary  parts  represent  the  frequencies  of  the  feedback 
controlled  vibration,  and  the  real  parts  represent  the  product  of  the 
daaping  ratio  and  the  natural  frequencies  of  the  undaaped  vibration. 
Notice  that  daaping  has  been  introduced  to  the  systea  through  the 
optinal  feedback  control,  as  evidenced  by  presence  of  nonzero  real  parts 
of  Eq  (150). 

The  nodal  natrix  of  this  systea  is  foraed  by  grouping  the 
eigenvectors  associated  with  the  eigenvalues  X  and  X  in  a 

d  ,  o  r ,  0 

natrix  fora  and 
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[♦]  * 


1.0000 

1.0000 

1.0000 

1.0000 

-0.9817 

+ 

0.13871  -0.9817 

- 

0. 1 387 i 

0.9890 

+ 

0.0986i 

0.9890 

- 

0.0986i 

0.3504 

- 

0.35661  0.3504 

+ 

0.35661 

0  4151 

- 

0. 4256i 

0.4151 

+ 

0. 4256i 

-0.2945 

♦ 

0.39871  -0.2945 

- 

0. 3987i 

0.4525 

- 

0.38001 

0.4525 

+ 

0. 3800 i 

-0.0044 

- 

0.2499i  -0.0044 

0. 2499i 

-0.0088 

- 

0. 3533i 

-0.0088 

+ 

0. 3533i 

0.0390 

+ 

0.24471  0.0390 

- 

0. 2447i 

0.0261 

- 

0. 3503i 

0.0261 

+ 

0. 3503i 

-0.0907 

- 

0.0860i  -0.0907 

+ 

0.0860i 

-0.1540 

- 

0. 1429 i 

-0.1540 

+ 

0. 1429i 

0.1099 

+ 

0.0718i  0.1009 

- 

0.0718i 

-0.1382 

- 

0. 15651 

-0.1382 

+ 

0. 1565i 

Fro*  the  social  aatrix,  one  say  construct  the  aodal  aatrices  associated 
with  the  traditional  state  vector  and  the  fractional  vector  of  the 
systea,  in  the  saae  Banner  as  indicated  in  the  first  part  of  this 


probleB.  For  the  current  feedback  control  scheae,  the  aatrices 
associated  with  the  state  vectors  are 


IV  s 

-0.0907  -  0.0860i 
0.1099  +  0.07181 
0.3504  -  0. 3566i 
-0.2945  ♦  0.39871 


-0.0907  +  0.0860i 
0.1009  -  0.0718i 
0.3504  +  0. 3566i 
-0.2945  -  0.3987i 


-0.1540  -  0.14291 
-0.1382  -  0.15651 
0  4151  -  0.42561 
0.4525  -  0.38001 


-0.1540  +  0.14291 
-0.1382  +  0. 15651 
0.4151  +  0.4256i 
0.4525  +  0. 3800i 


and 


[♦B]  - 
1.0000 

0.3504  -  0.35661 
-0.0044  -  0.24991 
-0.0907  -  0.08601 


1.0000 

0.3504  ♦  0. 3566i 
-0.0044  +  0.24991 
-0.0907  ♦  0.08601 


1.0000 

0  4151  -  0.42561 
-0.0088  -  0. 3533i 
-0.1540  -  0. 1 429 i 


1.0000 

0.4151  +  0.42561 
-0.0088  +  0. 3533i 
-0.1540  ♦  0. 1429i 
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The  state  transformation  aatrix  for  the  system  is 


<t>  *  [♦Al[*BJ't 


0  0 

1.2696  -3.2370 
0  1 

4.2675  -22.1595 


0 

4.2724 

0 

42.8966 


1 

0.1719 

0 

-37.4050 


This  state  transformation  matrix  permits  one  to  apply  the 
fractional  order  state  feedback  control  law  to  this  structure  and  the 
control  law  is  defined  as 


k  =  Ax  -  BG*z  (60) 

where  the  fractional  derivatives  of  the  system  response  at  mt  are  fed 
back  as  control  input  for  controlling  the  structure.  Defining  the 
fractional  state  gain  matrix,  G* ,  as 

G*  *  G*  (151) 

the  state  equations  for  fractional  order  state  feedback  control  can  be 


rewritten  as 


and 


-0.0371 

0.2337 

I -0.1198 

0.0034 

0 

0 

i  0 

0 

The  fractional  state  gain  aatrix,  G* ,  persits  the  direct  feedback  of 
fractional  states  for  the  control  of  the  structure.  The  feedback 
control  force  G*x  has  the  saae  control  effect  as  the  traditional  state 
feedback  control  force  G&. 


In  the  last  part  of  this  exaaple,  a  hypothetical  viscoelastic 
daaping  pad  is  introduced  at  »2  in  addition  to  the  control  force  at  st. 
The  systea  configuration  is  shown  in  Figure  5. 


Figure  5.  Viscoelastically-Daaped  2  D OF  Systea 


The  equations  of  aotion  for  the  systea  are  now 


[M]{w(t)>  +  [Kv]{A(t)}  ♦  [KE]{w(t)>  =  [D]{u(  t) } 


where 


[Kyi  = 


0 

0 


0 

0.1 


6  -2 

-2  6 


0  =  1/2 


(68) 


Before  one  aay  apply  the  optiaal  control  theory,  it  is  necessary  to 
develop  the  equivalent  daaping  and  equivalent  stiffness  aatrices  for  the 
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I 


systea  such  that  the  equations  of  notion  say  be  posed  in  the  integer 

i 

differential  fora  as 


[M]{w(t)>  ♦  rc*]<w(t)>  +  [K*l {w(t) }  =  [D] {u( t) }  (69) 

Following  the  procedure  identified  in  Section  3.2,  the  aodal 
acceleration,  velocity,  and  displaceaent  aatrices  for  the  hoaogeneous 
systea  aay  be  deterained  by  solving  the  fractional  order  equations  of 
aotion  of  the  systea 

D1/2[M]|d3/2  +  D1  ♦  D1/2  ♦  D° j { w ( t ) } 

-  [M][d3/2  ♦  D1  +  D1/2j{w(t)} 

+  [Kv]D1/2{w(t)}  +  [KE]D°{w(t)}  =  {0}  (153) 

through  the  use  of  the  solution  technique  of  the  expanded  equations  of 
notion.  The  expanded  equations  of  notion  for  Eq  (153)  are  now 
written  as 

D1/2[M]{w(tU  +  [K]{w(t) }  *  {0}  (154) 

where 


[0] 

[0] 

[0] 

[N] 

[0] 

[0] 

[N] 

[01 

(155) 

[0] 

in) 

[0] 

[01 

CM] 

[0] 

[0] 

[Kv]. 

67 


(156) 


’  [0] 

[0] 

-IM] 

[0] 

'N* 

[K]  = 

[0] 

-[M] 

[0] 

[0] 

-fMj 

[0] 

[0] 

[0] 

[°] 

[0J 

[0] 

[ke] 

and  (w(t)j  is  defined  by  Eq  (136).  The  modal  acceleration,  velocity, 


and  displacement  matrices:  [$#],  1^],  and  [♦  ],  are  then  determined 


from  the  eigenvalues  and  the  modal  matrix  of  the  system  posed  in  Eq 
(154)  and 


(4>J 


-1.4128  -  1.3219i  1.1894  +  1.19791 

1.4141  +  1.4206i  1.1888  ♦  1.12771 


[*  I 

V  J 


-0.3274  +  0,3531i  0.4174  -  0.42041 
0.3520  -  0.35361  0.3928  -  0.42001 


IV 


0.0883  +  0.081U  -0.1486  -  0.1454i 
-0.0884  -  0.0872i  -0.1484  -  0.13681 


Applying  the  above  matrices  to  Eq  (79)  and  (80),  the  equivalent  damping 
and  equivalent  stiffness  matrices  are  found  to  be 


[C*] 


0  0 
0.0034  0.0387 


[K*l 


6.0000  -2.0000 

-2.0112  6.1317 
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and  the  resonant  frequencies  of  the  viscoelastically  daaped  systea  are 

X*  S  -  0.0183  ±  4.0177i 

*'6  (157) 

X,  =  -  0.0204  ±  2. 8495i 

7,8 

By  posing  the  equations  of  notion  in  the  equivalent  fora,  Eq  (69). 
the  state  equations  can  now  be  developed  for  deteraining  the  optiaal 
gain  aatrix  for  the  feedback  control  of  this  viscoelastically  daaped 
structure,  following  the  saae  procedure  identified  in  Section  3.2.  The 
gain  aatrix  for  the  present  control  scheae  is 

G  »  [0,i02] 

’  I 

0.0066  -0.0048  j  0.1076  -0.0029 
*  0  0  I  0  0 

With  the  gain  aatrix  deterained,  the  optiaal  feedback  control  is 
then  introduced  into  the  systea  and  the  equations  of  notion  are  now 

[M]{w(tU  +  UC*]  +  [D][G2]){w(t)} 

+  UK*]  +  [D][Ot]){w(t)>  *  {0}  (84) 

For  a  *  1/2,  the  fractional  order  equations  of  notion  for  Eq  (84)  can  be 


written  as 


D1/2[M][d3/2  ♦  D1  +  D1/2  ♦  D°J {w(t ) > 

-  [Mj|d3/2  +  D1  +  D1/2J{w(t) } 

+  Dt/2([C*]  ♦  [D][Q21)[d1/2  ♦  D°] {w( t) } 

-  UC*]  +  [D][Q2])D1/2{w(t)} 

♦  UK*]  +  [D][G1])D°{w(t)]  =  {0}  (158) 

Eq  (158)  can  be  solved  using  the  solution  nethod  of  expanded  equations 
of  notion  described  in  Eq  ( 32 )-Eq  (36)  and  the  expanded  equations  of 
notion  for  Eq  (158)  are 

DUZ[M]{w(t)}  +  [K]{  w(  t ) }  =  (0)  (159) 

where 


[0] 

[0] 

[0] 

[Ml 

[M]  = 

[0] 

[0] 

[N] 

[0] 

(160) 

[0] 

[M] 

[0] 

[C* 

]+[D)[aj 

[Ml 

[0] 

[C*]+[0][02] 

[0] 

• 

’  [0] 

[0] 

- 

[M] 

[0] 

■ 

IK]  * 

[0] 

-[Ml 

[0] 

[0] 

(161) 

-[M] 

[0] 

-[C* 

1-[d][g2 

,] 

[0] 

[0] 

[0] 

[0] 

[K1 

r]  +  [D] [G  ] 

and  (w(t)}  is  defined  by  Eq  (136).  The  response  of  the  systen  can  be 
deternined  by  solving  the  expanded  equations  of  notion.  There  are  four 
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pairs  of  complex  conjugate  eigenvalues  to  Eq  (159) 


X4  ,  =  -  1.4047  ±  1. 4301i 

t.2 

X3  4  S  -  1.1781  ±  1.2093i 

Xc  =  1.4047  ±  1. 4301i 

s,6 

X  =  1.1781  ±  1.2093i 

7,8 

The  eigenvalues  X_  _  and  X  are  associated  with  the  aodal  behaviors  of 

5,6  7,8 

the  systea  and  the  aodal  frequencies  for  the  systea  are 


0.0718  ±  4.0177i 


X?  „  s  -  0.0745  ±  2.8495i 
7,8 


(162) 


Notice  that  there  is  an  increase  in  daaping  due  to  the  incorporation  of 
optiaal  feedback  control  with  passive  viscoelastic  daaping.  This  is 
evident  by  the  increase  in  aagnitude  for  the  real  parts  of  Eq  (162) 
coapare  with  the  real  parts  of  Eq  (157). 

The  aodal  aatrix  for  this  systea  is  foraulated  using  the  saae 
procedure  as  indicated  the  first  part  of  this  problea  and 


[♦]  * 


•0.9636  •  0.07261 

1.0000 

-0.3627  ♦  0.31761 
0.3496  -  0.35591 
-0.0138  +  0.240U 
-0.0044  -  0.24881 
0.0806  +  0.08881 
-0.0901  -  0.0854i 


-0.9636  ♦  0.07261 

1.0000 

-0.3627  -  0.3176i 
0.3496  *  0.35591 
-0.0138  -  0.240H 
-0.0044  ♦  0. 2488i 
0.0806  >  0.0888i 
-0.0901  +  0.0854i 


1.0000 

0.9688  +  0.043H 
0  4133  -  0. 4243i 
0.4187  -  0.39321 
-0.0092  -  0.35071 
0.0062  -  0.340U 
-0.1526  -  0 . 14111 
-0.1417  -  0.14321 


1.0000 

0.9688  -  0.0431i 
0.4133  +  0. 4243i 
0.4187  +  0. 3932i 
-0.0092  +  0. 3507i 
0.0062  t  0. 3401i 
-0.1526  +  0. 14111 
-0.1417  +  0. 1432i 
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The  nodal  Matrices  associated  with  the  traditional  state  vector  and  the 
fractional  state  vector  can  be  formulated  froa  the  systea  modal  matrix, 
following  the  procedure  described  in  the  first  part  of  this  problem  and 
they  are 


<V  = 

0.0806  +  0.0888i 
-0.0901  -  0. 0854i 
-0.3627  +  0.3176i 
0.3496  -  0. 3559i 


[♦.1  * 

-0.9636  -  0.07261 
-0.3627  +  0.3176i 
-0.0138  +  0.240U 
0.0806  +  0.08881 


0.0806  -  0.08881 
-0.0901  ♦  0.0854i 
-0.3627  -  0.3176i 
0.3496  +  0. 3559i 


-0.9636  +  0.07261 
-0.3627  -  0.31761 
-0.0138  -  0.24011 
0.0806  -  0.08881 


-0.1526  -  0. 14111 
-0.1417  -  0. 1432 i 
0  4133  -  0.4243i 
0.4187  -  0.3932i 


1.0000 

0  4133  -  0.42431 
-0.0092  -  0.35071 
-0.1526  -  0.14111 


-0.1526  +  0. 141 1 i 
-0.1417  +  0. 1432 i 
0.4133  +  0. 4243i 
0.4187  +  0.3932i 


1.0000 

0.4133  +  0.42431 
-0.0092  +  0.35071 
-0.1526  ♦  0.14111 


The  state  transformation  matrix  for  the  current  systea  is 


*  =  rojtV’1 


0  0  0  1 

1.2821  -3.2777  4.2951  0.1951 

0  10  0 

4.2822  -22.5219  43.7790  -38.3175 


This  transformation  matrix  may  now  be  introduced  into  the  state 
equations  to  incorporate  the  feedback  of  fractional  order  states,  and 
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the  new  state  equations  are 


k  *  Ax  -  BG*z 


(152) 


and 


-0.0186 

0.1887 

!  -0.1476 

0.1168 

0 

0 

i  o 

0 

The  foraulation  of  the  fractional  state  gain  aatrix  peraits  the 
direct  application  of  fractional  state  feedback  for  the  control  of 
structure  that  is  viscoelastically  daaped.  The  control  theory  reaaina 
unchanged  froa  the  undaaped  case  to  the  viscoelastically  daaped  case, 
while  the  coaposition  of  the  gain  aatrices  are  different  due  to  the 
variation  in  structural  coaponents. 
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Exaaple  Problem  2 


This  exaaple  is  part  of  an  ongoing  research  effort  here  at  AFIT  to 
demonstrate  the  application  of  fractional  order  state  feedback  control 
theory  in  the  active  control  of  large  flexible  space  structures.  The 
objective  of  this  exaaple  is  to  develop  the  optiaal  gain  aatrix  and  the 
state  transformation  aatrix  for  iapleaenting  fractional  order  state 
feedback  control  on  the  modified  Advanced  Beaa  Experiment  (ABE). 

The  modified  ABE  consists  of  a  six-foot  beaa,  with  the  attached  end 
welded  to  the  support  structure  and  a  circular  mounting  plate  attached 
to  the  free  end  of  the  beaa  (see  Figure  6).  The  physical  properties  of 
the  structure  is  listed  in  Table  1.  A  complete  description  of  the 
modified  ABE  setup  can  be  found  in  Reference  (25).  For  the  purpose  of 
this  analysis,  it  is  assumed  that  only  the  planar  bending  nodes  about 
the  Z-axis  would  be  present  in  any  response  and  no  other  bending  or 
torsional  modes  would  be  excited.  A  twenty  two  degree-of- freedom  finite 
element  model  of  the  inverted  cantilever  beaa  is  developed  for  this 
analysis.  The  circular  plate  is  modeled  as  a  point  mass  and  a  mass 
moment  of  inertia  combination  at  the  free  end  of  the  beaa.  The  attached 
end  connection  of  the  beaa  is  modeled  as  a  linear  spring  and  a 
rotational  spring  combination  to  better  represent  the  flexible  boundary 
condition.  The  two-spring  model  of  the  attached  end  also  allows  further 
fine  tuning  of  the  ABE  model  to  match  the  measured  baseline 
characteristics  of  the  beaa. 
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Table  1.  Structure  Physical  Properties 


Material  :  Aluainua 


Property  Description 

Value 

Unit 

Beaa  Length  (L) 

70.75 

in 

Y  Cross-Section  Width  (a) 

1.01 

in 

Z  Cross-Section  Width  (b) 

0.758 

in 

Cross-Section  Area  (A) 

0.7656 

,  2 
in 

Young’s  Modulus  (E) 

X 

00 

o 

H 

L06 

lbf/in2 

Shear  Modulus  (G) 

4.1  x  106 

lbf/in2 

Beaa  Density  (p) 

2.591  x 

10" 

4 

lbf -sec2/in 

Beaa  a ass  (aw ) 

D 

1.403  x 

10* 

2 

2 

lbf  sec  /in 

Y  Moaent  of  Inertia  (Iy) 

3.667  x 

10" 

2 

.  A 

in 

Z  Moaent  of  Inertia  ( I ^ ) 

6.508  x 

10" 

tm 

.  A 

in 

Torsional  Moaent  of  inertia  (K) 

7.913  x 

10" 

2 

.  A 

in 

Polar  Moaent  of  Inertia  (I  ) 

•  X 

1.865  x 

10" 

3 

2 

lbf -sec  /in 

Plate  Diaaeter  (d) 

12.0 

in 

Plate  Thickness  (t) 

1.0 

in 

Plate  Mass  (a  ) 

p 

2.847  x 

10" 

2 

2 

lbf -sec  /in 

Plate  X  Mass  Moaent  of  Inertia  (In  ) 

Dux 

0.5125 

lbf -sec2  -  in 

Plate  Y-Z  Mass  Moaent  of  Inertia  d0|||) 

0.2562 

lbf -sec2- in 
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«1  =  12.355  rad/sec  =  1.966  Hz 
<*>2  =  147.24  rad/sec  *  23.435  Hz 

and  they  are  in  cloae  agreement  with  the  aeaaured  resonances  of  the 
final  ABE  configuration  (25:41) 


=  1.97 

Hz 

1 , 

u 

=  23.43 

Hz 

2 , »••• 

In  this  exaaple,  it  is  desired  to  apply  the  fractional  order  state 
feedback  control  for  controlling  the  first  two  bending  nodes  of  the 
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beaa.  To  begin  the  analysis  of  this  problem,  the  bending  notion  of  the 
beaa  is  measured  with  acceleroaeters  at  positions  2,  5,  8,  and  11.  With 
traditional  state  feedback  control,  the  translational  displaceaents  and 
velocities  at  these  four  nodes  are  fed  back  to  the  controller  at 
position  5  to  provide  control  force  input  for  controlling  the  vibration 
of  the  beaa  (see  Figure  7).  Since  the  structural  control  of  the  beam 
is  based  on  the  nodal  response  at  the  four  nodal  positions,  it  is 
necessary  to  determine  the  optiaal  gain  aatrix  based  on  these  responses. 

The  equations  of  notion  for  the  systea  can  be  expressed  as 


[M] (w(t) }  +  [K] { w( t) }  =  [D] { u ( t ) } 


(165) 


where 


(w(t)}  *  { 


d 


0. 


'll 
011  > 


(166) 


[N]  and  [K]  are  the  aass  and  stiffness  aatrices  of  the  finite  eleaent 
aodel,  and  (w(t)}  is  a  displacement  vector  that  includes  both 
translational  and  angular  displaceaents  of  the  beaa  at  the  nodes.  The 
product  of  [D]  and  [u(t)J  represents  the  control  force  applied  to  the 
systea.  In  this  case  [0]  is  a  22  *  22  zero  aatrix  except  for  an  entry 
of  one  at  the  aatrix  position  D(9,9)  to  represent  the  control  force 
applied  at  nodal  position  5  of  the  beaa. 
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Position  9 


Figure  7.  ABE  Configuration  for  Traditional  State  Feedback 


The  analysis  of  the  gain  matrix  begins  with  the  formulation  of  the 


state  equations  of  the  system  with  optimal  control 

X  =  Ax  -  BGg  ( 1C7) 


where 


X  =  <  (168) 

Jw(t)}j 

X  is  the  state  vector  of  the  system  with  44  states,  that  includes  both 
displacement  and  velocity  at  all  degrees  of  freedom.  [A]  and  [B]  are 
constant  matrices  defined  by  Eq  (S3)  and  (54).  Q  is  the  optimal  gain 
matrix  for  the  present  system  with  full  state  feedback  control,  and  it 
is  derived  from  the  solution  of  the  Riccati  equation.  For  the  quadratic 
performance  index 


CD 

J  =  -jr  J ix(t)TQx(t)  +  uT(t)Bu(t)]dt  (169) 

o 


where 


Q 


[K]  [0] 

.10)  [M], 


(170) 


R  = 


1 

0 


0 

1 


(170) 
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the  optiaal  gain  aatrix,  6,  for  full  state  feedback  is  a  22  *  44  aatrix 
with  one  row  of  gain  elements  at  row  9  and  no  entries  at  other  rows. 

This  corresponds  to  the  full  state  feedback  control  setup  with  one 
control  force  applied  in  the  direction  of  the  translational  displaceaent 
at  node  5  of  the  beaa.  The  gain  elenents  at  row  9  of  the  full  state 
feedback  gain  aatrix  are 


G(9, 1:44)  = 
coluans  1  through  i 
[  0.3374  -2.7632 

coluans  10  through 
9.6753  1.1421 

coluans  19  through 
-3.9078  -2.4949 
coluans  28  through 
-0.0087  -0.0028 

coluans  37  through 
0.0298  0.0058 

coluans  41  through 
-0.0086  0.0255 


-0.6534  0.5310 

16 

-2.4241  8.3482 

24 

0.2251  0.4790 

32 

-0.0096  0.0011 

40 

0.0119  0.0075 

44 

0.2232  0.1307] 


-0.9246  -0.4154 

-0.4293  4.1572 

1.3739  -2.6400 

-0.0103  -0.0111 

0.0189  0.0367 


-3.9616  -8.0019 

0.9121  7.4925 

0.0122  -0.0046 

0.1595  0.0207 

-0.0413  -0.0406 


Foraulating  the  gain  aatrix  G  *  j G2 ]  and  apply  the  feedback  control 
to  the  systen,  the  equations  of  notion  of  the  beaa  nay  now  be  written  as 


[Mj{w(tU  +  [D][G2]{w(t)}  +  UK]  +  [D][Gtl){w(t)}  =  {0}  (56) 
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The  systea  response  of  the  full  state  feedback  controlled  beaa  can  be 
deterained  by  solving  Eq  (56)  in  first  order  fora 


d’emHwUH  +  [K] [w( t) }  =  (0) 


(172) 


where 


[MJ  = 

[0] 

(Ml 

(173) 

.  (Ml 

(dkg23  . 

[K]  = 

-[M] 

(01 

(174) 

.  col 

[Kl  +  IDHG  3  . 

(w(t)} 


a  < 


Dl{w(t)) 

(w(t)) 


} 


(175) 


and  the  frequencies  for  the  first  four  Bodes  of  vibration  are 


=  -  0.9725  ±  12.355i 
w  a  -  7.1178  *  147. 24i 

(176) 

«3  *  -  5.0212  ±  381.561 
a  -  3.4060  *  638.531 

The  daaping  effect  introduced  by  the  feedback  control  are  evident  as 
shown  by  the  negative  real  parts  of  the  systea  responses. 

Since,  for  this  analysis,  the  notion  of  the  beaa  is  only 
represented  by  the  translational  displaceaents  and  velocities  at  the 
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nodal  positions  2,  5,  8,  and  11,  the  application  of  full  state  feedback 
control  using  the  complete  44  states  of  the  systea  is  not  very 
practical.  It  is  desirable  to  develop  an  equivalent  control  force  that 
has  the  saae  effect  as  the  full  state  feedback  control  while  using  a 
reduced  set  of  states,  xR,  and  xR  is  defined  in  this  case  as 


It  can  be  shown  that  there  exists  a  relation  between  the  full  state 
vector,  g,  and  the  reduced-state  vector,  jR,  through  a  state  reduction 
transforaation  aatrix,  [$R]. 

Recalling  froa  structural  aodal  analysis  that  the  systea  response 
aay  be  expressed  as  a  product  of  the  aodal  aatrix  and  the  aodal 
coordinate  vector.  Express  the  full  state  vector,  x,  and  the  reduced 
state  vector,  xR,  is  terns  of  the  their  respective  nodal  aatrix  and  the 
aodal  coordinates 

X  a  [♦'Hr?)  (178) 

XR  *  If] in)  (179) 

where  [♦'  ]  and  [♦"!  are  the  associated  aodal  aatrices  for  &  and  x„> 


respectively.  The  row  elements  of  [$' ]  and  [$"]  are  those  row  elements 
of  the  full  aodal  aatrix,  [♦] ,  of  the  full  state  feedback  control  system 
in  Eq  (172)  associated  with  the  states  in  x  and  xR*  The  full  aodal 
aatrix,  [♦] ,  is  a  square  aatrix  consist  of  conjugate  eigenvectors  of  the 
systea  in  Eq  (172).  It  is  desired  for  the  equivalent  control  force  to 
produce  the  saae  control  effect  over  the  first  four  nodes  of  vibration 
as  the  full  state  feedback  control  force,  thus  only  those  eigenvectors 
of  the  full  aodal  aatrix,  [♦] ,  that  are  associated  with  the  first  four 
nodes  are  used  in  the  foraulaticn  of  [♦' ]  and  [♦"],  where  [0' ]  is  a 
44  x  8  aatrix  and  [4"j  is  a  8  *  8  aatrix.  Froa  Eq  (179)  one  aay  define 
(H)  as 

{n>  *  (180) 

Substitute  Eq  (180)  into  Eq  (178)  yields 

X  ■[♦'It*'  f'xR  (181) 

Defining  the  state  reduction  transfornation  aatrix,  (^R]»  as 

[*R]  *  [♦'! [♦"!''  (182) 

and  Eq  (181)  aay  be  posed  as 

X  *  (183) 
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and  Eq  (183)  illustrates  the  relation  between  the  full  state  vector  and 
the  reduced  state  vector  for  the  first  four  aodes  of  vibration. 

Substitute  Eq  (183)  into  Eq(167)  and  the  state  equations  can  be 
written  as 

k  =  Ax  -  BG[0r]£r  (184) 

Define  the  reduced  gain  aatrix  as 

Gr  =  G[*r]  (185) 

and  Eq  (184)  can  be  written  as 

X  =  Ax  -  BGr*r  (186) 

where  GR  is  a  22  x  8  matrix  with  one  row  of  gain  elements  at  row  9.  Eq 
(186)  is  the  state  equations  with  an  equivalent  control  force  being 
applied  in  the  translational  direction  at  nodal  position  5,  and  the 
equivalent  control  force  is  developed  based  on  the  states  of  the  reduced 
state  vector  defined  in  Eq  (177).  The  gain  elesents  at  row  9  of  the 
reduced  gain  aatrix  of  Eq  (185)  are 

GR(9,1:8)  * 
columns  1  through  8 

(-1.1477  2.4025  -1.6402  0.4945  -0.0482  0.2196  -0.0371  0.2224] 
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and  the  gain  elements  represent  the  gains  to  be  applied  to  the  states  of 
the  reduced  state  vector  in  Eq  (177)  for  generating  an  equivalent 
control  force  for  controlling  the  first  four  aodes  of  vibration. 

To  introduce  the  gains  for  the  reduced  state  vector  into  the  state 
equations,  Eq  (167),  one  simply  expands  GR  columnwise  by  incorporating 
zero  coluan  vectors  into  those  columns  that  correspond  to  the  states  of 
the  full  state  vector  not  considered  for  the  feedback.  This  expansion 
process  results  in  a  gain  matrix,  GR£,  with  the  saae  gain  eleaents  froa 
the  reduced  gain  matrix  expanded  across  row  9  of  the  matrix  GR£.  The 
gain  elements  at  row  9  of  the  matrix  GR£,  are 

<>,.(9.1:44)  = 

columns  1  through  8 


[  0.0 

0.0 

-1.1477 

0.0 

0.0 

0.0 

0.0 

0.0 

columns  9  through  16 

2.4025 

0.0 

0.0 

0.0 

0.0 

0.0 

-1.6402 

0.0 

columns  17 

through 

24 

0.0 

0.0 

0.0 

0.0 

0.4945 

0.0 

0.0 

0.0 

columns  25 

through 

32 

-0.0482 

0.0 

0.0 

0.0 

0.0 

0.0 

0.2196 

0.0 

columns  33 

through 

40 

0.0 

0.0 

0.0 

0.0 

-0.0371 

0.0 

0.0 

0.0 

columns  41 

through 

44 

0.0 

0.0 

0.2224 

0.0  ] 
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The  state  equations  for  the  systea  with  reduced  state  feedback  is 


k  =  Ax  -  bgrex 


(187) 


and  the  product  of  GREx  has  the  saae  control  effect  on  the  first  four 


■odes  as  G„x  .  Formulate  the  new  gain  aatrix  as  G__  =  [G. 


J*  the 


RR  -  -  -  o  RE  "~RE1 ' “RE2 ' 

equations  of  notion  for  the  systea  with  reduced  feedback  control  can  be 
written  as 

tM]{w(t)}  ♦  [DHGRE2]{w(t)}  ♦  ([K]  +  tD] CGRE11  ){w(t) }  =  {0}  (188) 

The  systea  response  aay  be  deterained  by  solving  the  expanded  equations 
of  notion  using  the  procedure  identified  in  Eq  (172)-(175),  and  the 
frequencies  for  the  first  four  nodes  of  vibration  are 


*  -  0.9725  ±  12. 35 5 i 
u  =  -  7.1178  ±  147. 24i 

*  (189) 

w3  s  -  5.0212  ±  381.561 

=  -  3.4060  ±  638.531 

i 

Notice  that  the  frequencies  for  the  first  four  nodes  in  Eq  (189),  using 
the  reduced  state  feedback  control,  are  identical  to  the  response  of  the 
^  full  state  feedback  control  in  Eq  (176). 

At  this  point  of  the  analysis,  it  is  desired  to  introduce 
fractional  order  state  feedback  control  in  place  of  the  reduced  state 
feedback  control.  The  acceleroaeters  at  positions  2,  5,  8,  and  11  are 


8 


replaced  by  two  accelerometers  at  positions  5  and  8,  plus  the  associated 
electronic  fractional  order  differentiator  circuits  of  order  a  =  1/2  . 
The  fractional  derivatives  of  the  translational  displacements  at  these 
two  nodes,  together  with  their  displacements  and  velocities,  are  fed 
back  to  the  controller  to  provide  control  force  input  for  controlling 
the  structure  (see  Figure  8).  The  application  of  the  fractional  order 
state  feedback  control  theory  in  this  example  will  reduce  the  sensor 
requirement  from  four  to  two,  while  still  maintaining  the  same  control 
authority  over  the  system  as  the  traditional  state  feedback  control. 

To  introduce  fractional  order  state  feedback  control  to  the  beam, 
one  needs  to  solve  the  fractional  order  equations  of  motion  of  the 
system  in  the  expanded  format  and  formulate  the  state  transformation 
matrix  following  the  procedure  identified  in  Eq  (32)-Eq  (48).  For 
a  =  1/2,  the  fractional  order  equations  of  motion  of  Eq  (188)  are 

D1/2[M][d3/2  +  Dl  +  D1/Z  +  D° j { w ( t ) } 

-  [M][d3/2  +  D1  ♦  D1/2]{w(t)} 

♦  D1/2[D][0RE2][d1/2  ♦  D° j { w ( t )} 

'  CDlCGRE2]I>,/2fw(t:)> 

+  ( [K]  ♦  [D][GRE1J)D°{w(t)}  =  {0}  (190) 

and  Eq  (190)  can  be  solved  using  the  solution  method  of  the  expanded 
equations  of  motion  described  in  Eq  (32)-Eq  (36).  The  expanded 
equations  of  motion  for  Eq  (190)  are 

D1/2[M]{w(t) )  +  (K]{w(t) }  *  {0}  (191) 
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Figure  8.  ABE  Configuration  for  Fractional  State  Feedback 
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1 


(193) 


(194) 


and  (w(t)}  for  this  systea  is  defined  in  Eq  (166).  The  nodal  satrix  of 
the  systea,  [4] ,  is  composed  of  the  coaplex  conjugate  pairs  of 
eigenvectors  associated  with  the  conjugate  eigenvalues  on  the  principal 
sheet  of  the  Bieaann  surface.  Since  only  the  first  four  nodes  of 
vibration  are  controlled  by  the  equivalent  control  force,  only  the 
eigenvectors  correspond  to  the  first  four  nodes  of  vibration  will  be 
used  in  foraulating  [♦]  and  the  resulting  aatrix  is  88  *  8.  It  is 
desired  to  deternine  the  state  transforaation  aatrix,  <P,  for  the 
application  of  feeding  back  the  fractional  states  of  responses  at  nodal 
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'i 


positions  5  and  8.  Invoke  the  expansion  theorem  and  the  full  state 
vector  aay  be  written  as 


x  =  [*A]{n> 


for 


(195) 


.  (w(t)}l 


(196) 


[0  3  -  a  44  x  8  aatrix  consist  of  row  vectors  of  [♦]  associated 
with  the  states  of  the  full  state  vector  x 


The  fractional  state  vector  aay  be  written  in  the  fora 


X  =  [«>B]{n> 


for 


r  'i 

d3/2< 

y5l 

► 

•V 

D  1  « 

Js  ' 

► 

r«. 

-.1/2 

y 

D 

*5 
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ys 

► 

k  y8  , 

(197) 


(198) 


[♦ ]  a  a  8  k  8  aatrix  consist  of  row  vectors  of  [♦}  associated 

B 

with  the  states  of  the  fractional  state  vector  y 
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Eq  (203)  is  the  state  equations  for  the  feedback  of  reduced  fractional 

order  state  control,  where  G*g  is  a  22  *  8  satrix  with  one  row  of  gain 

«•» 

eleaents  at  row  9.  The  gain  elesent  at  row  9  of  GRE  are 


92 


G*e(9,1:8)  = 
coluans  1  through  8 

[-0.0226  0.0177  1.1198  -0.5565  -15.9949  9.6436  100.1118  -49.1084] 

This  gain  aatrix  persits  the  fractional  states  froa  positions  5  and  8  to 
be  fed  back  for  the  controlling  of  the  first  four  vibrational  Bodes  of 
the  structure. 

Theoi otically,  the  the  control  effect  provided  the  product  of  the 
gain  aatrix  G*R£  and  the  fractional  states  observed  at  positions  5  and  8 
for  controlling  the  first  four  nodes  of  vibration  of  the  beaa  is 
equivalent  to  the  control  effect  provided  by  the  product  of  the  gain 
aatrix  GRE  and  the  traditional  states  observed  at  positions  2,  5,  8,  and 
11.  But  physically  only  the  first  two  nodes  of  vibration  of  the  beaa 
can  be  accurately  controlled  by  this  control  systea.  This  is  due  to  the 
liaitation  of  the  finite  eleaent  aodel  of  the  beaa,  of  which  only  the 
first  two  nodes  of  vibration  can  be  aodeled  with  relative  accuracy. 

The  next  step  in  this  experiaent  is  to  incorporate  the  gain  aatrix 
to  the  fractional  order  state  feedback  control  of  the  beaa  and 
physically  deaonstrate  the  validity  of  the  control  theory.  This 
validation  process  is  currently  underway  at  AFIT  and  the  results  of  this 
experiaent  will  becoae  available  in  the  near  future. 
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VI.  Conclusion 


The  fractional  order  state  feedback  control  theory  is  developed  by 
incorporating  the  state  transformation  matrix  into  the  traditional 
state-space  equations  of  the  structure  under  control.  The  state 
transformation  matrix  between  the  traditional  states  and  the  fractional 
states  can  be  formulated  by  establishing  the  associated  modal  matrices 
from  the  eigenstructures  of  the  expanded  equations  of  motion  of  the 
system.  The  optimal  control  gains  are  determined  by  solving  the  Riccati 
equation  for  the  linear  quadratic  regular  problem.  The  gains  are  then 
combined  with  the  state  transformation  matrix  to  form  the  necessary 
fractional  state  gains  for  the  fractional  order  state  feedback  control 
for  the  system. 

The  use  of  fractional  order  state  feedback  control  appears  to  hold 
value  in  providing  an  alternative  form  of  an  observer  which  permits  one 
to  use  fewer  sensors  at  limited  locations,  while  still  be  able  to 
provide  the  equivalent  observability  and  control  authority  over  the 
structure  as  the  traditional  state  feedback  control  scheme. 

The  next  step  in  the  development  of  fractional  order  state  feedback 
control  is  to  test  the  theory  by  experiment.  The  current  research 
effort  on  the  modified  ABE  and  the  development  of  a  prototype  sensor  for 
fractional  derivative  motion  should  provide  insight  on  the  application 
of  the  theory  in  real-time  control  situation.  Additional  research 
efforts  are  still  needed  in  determining  the  optimal  number  of  sensors 
required,  the  optimal  sensor/controller  locations  on  the  structure  being 
controlled,  and  the  optimal  fractional  states  for  feedback  structural 
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control  in  order  to  further  expand  the  application  of  the  theory. 
Furthermore,  advanced  electronic  fractional  order  differentiator 
circuitries  are  needed  to  advance  the  development  and  performance  of  the 
fractional  order  state  control  system. 
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